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Abstract: The computation of exclusive QCD jet observables at higher orders requires 
a method for the subtraction of infrared singular configurations arising from multiple ra- 
diation of real partons. We present a subtraction scheme relevant for NNLO perturbative 
calculations in e + e~ — * jets. The building blocks of the scheme are antenna functions 
derived from the matrix elements for tree-level 1 — ► 3 and 1 — > 4 and one-loop 1 — ► 3 
processes. By construction, these building blocks have the correct infrared behaviour when 
one or two particles are unresolved. At the same time, their integral over the antenna 
phase space is straightforward. As an example of how to use the scheme we compute the 
NNLO contributions to the subleading colour QED-like contribution to e + e~ — * 3 jets. To 
illustrate the application of NNLO antenna subtraction for different colour structures, we 
construct the integrated forms of the subtraction terms needed for the five-parton and four- 
parton contributions to e + e~ — * 3 jets at NNLO in all colour factors, and show that their 
infrared poles cancel analytically with the infrared poles of the two-loop virtual correction 
to this observable. 
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1. Introduction 

Experimental measurements of jet production observables are among the most sensitive 
tests of the theory of Quantum Chromodynamics (QCD), and yield very accurate determi- 
nations of QCD parameters [1], especially of the strong coupling constant a s . At present, 
the precision of many of these determinations is limited not by the quality of the experi- 
mental data, but by the error on the theoretical (next-to-leading order, NLO) calculations 
used for the extraction of the QCD parameters. To improve upon this situation, an exten- 
sion of the theoretical calculations to next-to-next-to-leading order (NNLO) is therefore 
mandatory [2]. 

In the recent past, many ingredients to NNLO calculations of collider observables have 
been derived, including the universal three-loop QCD splitting functions [3] which govern 
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the evolution of parton distribution functions at NNLO. The massless two- loop 2^2 and 
1 — > 3 matrix elements relevant to NNLO jet production have been computed [4, 5] using 
several innovative methods [6], and are now available for many processes of phenomeno- 
logical relevance. The one-loop corrections to 2 — > 3 and 1 — > 4 matrix elements have been 
known for longer [7,8] and form part of NLO calculations of the respective multi-jet observ- 
ables [9,10]. These NLO matrix elements naturally contribute to NNLO jet observables of 
lower multiplicity if one of the partons involved becomes unresolved (soft or collinear) [11]. 
In these cases, the infrared singular parts of the matrix elements need to be extracted and 
integrated over the phase space appropriate to the unresolved configuration to make the 
infrared pole structure explicit. Methods for the extraction of soft and collinear limits 
of one-loop matrix elements are worked out in detail in the literature [11-16]. As a final 
ingredient, the tree level 2^4 and 1 — > 5 processes also contribute to (2 — > 2)- and 
(1 — ► 3)-type jet observables at NNLO. These contain double real radiation singularities 
corresponding to two partons becoming simultaneously soft and/or collinear [17-20]. To 
determine the contribution to NNLO jet observables from these configurations, one has 
to find two-parton subtraction terms which coincide with the full matrix element and are 
still sufficiently simple to be integrated analytically in order to cancel their infrared pole 
structure with the two- loop virtual and the one- loop single- unresolved contributions. In 
the past, such configurations were only dealt with on a case- by-case basis in the context 
of specific calculations [17,21-24], while no general method was available. Several meth- 
ods have been proposed recently to accomplish this task [25-29]. Up to now, only one 
method has been fully worked through for observables of physical interest:the sector de- 
composition algorithm [30,31]. In this method, both phase space and loop integrals are 
analytically decomposed into their Laurent expansion in dimensional regularisation, and 
the coefficients of the expansion are numerically integrated. Results have been obtained 
for e + e~ — > 2j [32], pp —>■ H + X [33] and most recently for muon decay [34] at NNLO. In 
contrast to all other approaches, in the sector decomposition method one does not have to 
integrate the subtraction term analytically. 

Following a number of pioneering calculations [35-37], infrared subtraction of real 
radiation singularities is well understood at NLO, where several generic process independent 
methods exist [38,40,41]. One of these methods is the so-called antenna subtraction [10,41], 
which derives NLO subtraction terms from three-parton tree-level matrix elements which 
naturally encapsulate all singular limits due to unresolved emission of a single parton 
between two colour-connected [42,43] hard partons (radiators). In [44], we described the 
construction of NNLO subtraction terms for e + e~ — > 2j based on full four-parton tree- 
level and three-parton one-loop matrix elements, which can be integrated analytically over 
the appropriate phase spaces [31]. If normalised appropriately, these full four-parton tree- 
level and three-parton one-loop matrix elements can be interpreted as antenna functions at 
NNLO. These NNLO antenna subtraction terms derived from four-parton matrix elements 
with a hard quark-antiquark pair in [44] were used subsequently [45] to compute the over- 
correction to e + e~ — > 3j at NNLO. 

The calculation presented in [44] included only antenna functions involving radiation 
off a hard quark-antiquark pair. We derived colour-ordered matrix elements (and thus 
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NNLO antenna functions) for the two other partonic configurations (hard quark-gluon ra- 
diators and hard gluon-gluon radiators) by using appropriate effective Lagrangian densities 
coupling an external current to a gluino-gluon [46] or gluon-gluon [47] system. 

In the present paper, we describe the antenna subtraction method at NNLO in detail 
and provide an algorithm for the construction of antenna subtraction terms at NNLO 
from the basic three-parton and four-parton antenna functions. For the sake of clarity, 
we restrict ourselves to the kinematical situation of a colour-neutral particle decaying into 
coloured final state partons, as relevant to e + e~ — ► m jets at NNLO. The basic structure of 
the NNLO antenna subtraction can however be carried over to configurations with partons 
in the initial state. 

The paper is structured as follows. In Section 2 we describe the antenna subtraction 
method at NLO and NNLO, and explain the construction of the subtraction functions 
for tree-level and one-loop real radiation matrix elements up to NNLO accuracy from 
the basic three- and four-parton antenna functions. All antenna functions are derived 
and integrated over their appropriate final state phase spaces in the following sections. 
Section 3 establishes the notation, and lists the different possible antenna configurations. 
In Section 4, we introduce the colour-ordered infrared singularity operators [48] appearing 
in the integrated NLO and NNLO antenna functions. The quark-antiquark, quark-gluon 
and gluon-gluon antenna functions in both unintegrated and integrated forms are then 
derived in Sections 5-7. Their behaviour in all single and double unresolved limits is 
summarised in Section 8. As a first non-trivial application of our method, we document 
the calculation of the subleading colour NNLO QCD correction to the e + e~ — > 3j cross 
section (which was already reported briefly in [45]) in Section 9. To illustrate the generality 
of the method, and to outline future applications of it, we reconstruct the infrared pole 
structure of the integrated NNLO subtraction terms for e + e~ — > 3j in all colour factors, 
and show that these cancel the explicit infrared poles of the two-loop virtual corrections to 
this observable in Section 10. Finally, Section 11 contains our conclusions and an outlook 
on applications and extensions of the method presented here. 

2. Infrared subtraction terms 

To obtain the perturbative corrections to a jet observable at a given order, all partonic 
multiplicity channels contributing to that order have to be summed. In general, each par- 
tonic channel contains both ultraviolet and infrared (soft and collinear) singularities. The 
ultraviolet poles are removed by renormalisation, however the soft and collinear infrared 
poles cancel among each other when all partonic channels are summed over [49]. 

While infrared singularities from purely virtual corrections are obtained immediately 
after integration over the loop momenta, their extraction is more involved for real emission 
(or mixed real-virtual) contributions. Here, the infrared singularities only become explicit 
after integrating the real radiation matrix elements over the phase space appropriate to the 
jet observable under consideration. In general, this integration involves the (often iterative) 
definition of the jet observable, such that an analytic integration is not feasible (and also 
not appropriate). Instead, one would like to have a flexible method that can be easily 
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Figure 1: Illustration of colour connection. 



adapted to different jet observables or jet definitions. Therefore, the infrared singularities 
of the real radiation contributions should be extracted using infrared subtraction terms. 
The crucial points that all subtraction terms must satisfy are that (a) they approximate 
the full real radiation matrix elements in all singular limits and (b) are still sufficiently 
simple to be integrated analytically over a section of phase space that encompasses all 
regions corresponding to singular configurations. Note that the subtraction terms should 
also be local and should fully account for the limit they are aimed at without introducing 
spurious infrared singularities in other limits. 

In this paper, we shall restrict ourselves to the kinematical situation of a massive, 
colour-neutral particle decaying into massless coloured partons, such as multijet production 
in electron-position annihilation: 

e + e~ — ► m jets . 

In this situation, infrared singularities appear only due to final state radiation. In the more 
general partonic scattering kinematics, coloured partons are also present in the initial state, 
yielding initial state and fixed initial/final state singularities. The method presented here 
can, like the antenna factorisation method at NLO [41], be extended to these situations. 
This extension is beyond the scope of this paper, and will only be commented on briefly in 
Section 11 below. 

To specify the notation, we define the tree-level n-parton contribution to the m-jet 
cross section (for tree level cross sections n = m; we leave n ^ m for later reference) in d 
dimensions by, 

da B = M £d*„(pi, ...,Pn, q)^~ \M n ( Pl , . . . , Pn )\ 2 Jg>(pi, ■ ■ ■ ,Pn). (2.1) 

the normalisation factor Af includes all QCD-independent factors as well as the dependence 
on the renormalised QCD coupling constant a s , J2 n denotes the sum over all configurations 
with n partons, d<E> n is the phase space for an n-parton final state with total four-momentum 
qV- in d = 4 — 2e space-time dimensions, 

while S n is a symmetry factor for identical partons in the final state. |.M n | 2 denotes 
a squared, colour-ordered tree-level n-parton matrix element, where particle 1 is colour 
connected to particle 2 which is colour connected to particle 3 and so on as illustrated in 
Figure 1. Contributions to the squared matrix element which are sub leading in the number 
of colours can equally be treated in the same context, noting that these subleading terms 
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yield configurations where a certain number of essentially non-interacting particles are 
emitted between a pair of hard radiators. By carrying out the colour algebra, it becomes 
evident that non-ordered gluon emission inside a colour-ordered system is equivalent to 
photon emission off the outside legs of the system [18,42]. For simplicity, these subleading 
colour contributions are also denoted as squared matrix elements |.M m | 2 , although they 
often correspond purely to interference terms between different amplitudes. 

The precise definition depends on the number and types of particles involved in the 
process. However, all colour orderings are summed over in J2m w ith the appropriate colour 

(n) 

weighting. The jet function J m defines the procedure for building m jets out of n partons. 

(n) 

The main property of J m is that the jet observable defined above is collinear and infrared 
safe as explained in [39,40]. In general jffl contains 6 and S- functions. jffl can also 
represent the definition of the n-parton contribution to an event shape observable related 
to m-jet final states. 

From (2.1), one obtains the leading order approximation to the m-jet cross section by 
integration over the appropriate phase space. 

da LO = [ da B . (2.3) 

Depending on the jet function used, this cross section can still be differential in certain 
kinematical quantities. 



2.1 NLO infrared subtraction terms 

At NLO, we consider the following m-jet cross section, 



da nlo 



(do-§LO 



dtTNLo) + 



/ 



da 



V 

NLO 



(2.4) 



The cross section da NLO has the same expression as the Born cross section da NLO (2.1) 
above except that m — > m + 1, while da NLO is the one- loop virtual correction to the m- 
parton Born cross section da B . The cross section da NLO is a (preferably local) counter-term 
for d<j B L0 . It has the same unintegrated singular behaviour as da NLO in all appropriate 
limits. Their difference is free of divergences and can be integrated over the (m + l)-parton 
phase space numerically. The subtraction term da NLO has to be integrated analytically 
over all singular regions of the (m + l)-parton phase space. The resulting cross section 
added to the virtual contribution yields an infrared finite result. 

A systematic procedure for finding NLO infrared subtraction terms is the antenna 
formalism introduced in [10,41]. The antenna subtraction terms are obtained as sum of 
antennae: 

daNLO = N Y] d^m+l (Pi , • • • , Pm+l ! <?) TT^— 

-. <->m+l 
m+1 

X Yl X ijk \M m (pi, . . . ,PI,PK, ■ ■ ■ ,Pm+l)| 2 4r } (Pl, • • • ,PI,PK, • • • ,Pm+l) , (2-5) 
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Figure 2: Colour connection of the partons showing the parent and daughter partons for the single 
unresolved antenna. 

such that, 

&°NLO ~ dcr NLO = N d$ m+l(Pl, • • • , Pm+l', <l) 

\M m +l(pi, ■ ■ ■ ,Pm+l)\ 2 J^ +1) (PU ■ ■ ■ ,Pm+l) 
~ ^2 X? jk \M m (pi, . . . ,PI,PK, ■ ■ ■ ,Pm+l)\ 2 J^{Pl, • • • ,PI,PK, • • • ,Pm+l) 

(2.6) 



The subtraction term involves the m-parton amplitude depending only on the rede- 
fined on-shell momenta pi, ■ ■ ■ ,Pi,Pk, ■ ■ ■ ,Pm+i where pi,pK are linear combinations of 
PiiPjiPk while the tree antenna function X^- k depends only on Pi,Pj,Pk- X^- k describes all 
of the configurations (for this colour-ordered amplitude) where parton j is unresolved. This 
occurs because of the particular factorisation properties of colour-ordered amplitudes. In 
particular, when particle j is a gluon and becomes soft, the colour-ordered matrix element 
undergoes QED-like factorisation so that, 

\M m+ l(pi, . • • ,Pi,Pj,Pk, • • • ,Pm+l)| 2 -> —\M m (pi, ...,pi,PK,--- ,p m +l)\ 2 - 

In this limit pi = pi and px = Pk- Together particles / and K form a colour connected 
hard antenna that radiates particle j. In doing so, the momenta of the radiators change to 
form particles i and k. The type of particle may also change. For example when particle 
j is a gluon, the colour factor is modified at the amplitude level by, 

as illustrated in Figure 2. The antenna factorisation of squared matrix element and phase 
space can be illustrated pictorially, as displayed in Figure 3. 

The antenna approach described here is closely related to the commonly used dipole 
factorisation formalism derived by Catani and Seymour [40]. Here, the corresponding 
subtraction term for the same colour-ordered squared matrix element is given by 



&°nlo = N ^2 d$ m+ i(pi, . . . ,p m +i;q)-^- 

— 



m+l 



m+1 



^2 V ij,k \M m {(pi, . . . ,Pij,Pk, • • • ,Prn+l)\ 2 J^HPI, ■ ■ ■ ,Pij,Pk, • • • ,Pm+l) 
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Figure 3: Illustration of NLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + l)-particle phase space. The term in square brackets repre- 
sents both the antenna function X^ k and the antenna phase space d$x ijk - 

+ ^2 V kj,i \M m ((pi, ■ ■ ■ ,Pi,Pkj, • • • ,Pm+l)\ 2 J^(Pl, ■ ■ ■ ,Pi,Pkj, • • • ,Pm+l) ■ 

(2.7) 

In the first term, the dipole contribution involves the m-parton amplitude which only 
depends on the redefined on-shell momenta p\, . . . ,Pij,Pk, ■ ■ ■ ,Pm+i and the dipole function 
Djj^ which depends on pi^p^p^. The momenta pi, pj and p^ are respectively the emitter, 
unresolved parton and the spectator momenta corresponding to a single dipole term. In 
the second term, the role of emitter and spectator are exchanged. The redefined on- 
shell momenta pij,pk (fikj,Pi) are different linear combinations of pi, pj and p\ for each 
dipole. In the antenna approach, the momentum mapping would be the same for each 
dipole contribution and the two terms combine to form the tree antenna, X^ k . The two 
dipoles combining to an antenna have a common unresolved parton, and contain the two 
possible emitter/spectator combinations. In the antenna language, emitter and spectator 
act as radiators. Note that we can always choose to divide the antenna and use different 
momentum maps for the two parts. 

The jet function J^f 1 in (2.6) does not depend on the individual momenta Pi, pj and 
Pk, but only on pi,pk- One can therefore carry out the integration over the unresolved 
dipole phase space appropriate to pi , pj and pk analytically, exploiting the factorisation of 
the phase space, 

d$ m+ i(pi, . . . ,p m+ i;q) = d$ m (pi, . . . ,Pi,Pk, ■ ■ ■ ,Pm+i;q) ■ d$x ijk (Pi,Pj,Pk;Pi +Pk) ■ 

(2.8) 

The NLO antenna phase space d$x ijk is proportional to the three-particle phase space, 
as can be seen by using m = 2 in the above formula and exploiting the fact that the 
two-particle phase space is a constant, 

P 2 = J dd> 2 = 2- 3+2 ^- 1+e ^f|y ( q 2 r ■ (2-9) 

such that 

d$ 3 = P2 d$x ljk ■ (2.10) 



-8- 



For the analytic integration, we can use (2.8) to rewrite each of the subtraction terms 
in the form, 



\M m \ 2 Jt ] d^mj <KI'y , A') 



where \M m \ 2 , jffl and d$ m depend only on p 1 , , . . . ,pi,Pk, ■ ■ ■ ,Pm+i and d$x ijk and Xf jk 
depend only on Pi,Pj,Pk- The analytic integral of the subtraction term is therefore defined 
as the antenna function integrated over the fully inclusive antenna phase space, normalised 
appropriately, 

Xijk^ijk) = (Svr 2 (4vrr ^) J d$ Xijk X% k . (2.11) 

This integration is performed analytically in d dimensions to make the infrared singu- 
larities explicit and added directly to the one-loop m-particle contributions. The factor 
(8-7T 2 (47r) _e e n ) in the above equation is related to the normalisation of the renormalised 
coupling constant, and its relation to the bare coupling parameter g = ^/47tqo appearing 
in the QCD Lagrangian density: 



where 



(2.12) 



S e = (47r) e e _e7 with Euler constant 7 = 0.5772 . . . 

and p,Q is the mass parameter introduced in dimensional regularisation to maintain a di- 
mensionless coupling in the bare QCD Lagrangian density; (5q and (5\ are the first two 
coefficients of the QCD /3-function: 

UN - 2N F 34iV 3 - 13N 2 N F + 3N F 
Po = g , Pi = m , (2.13) 

with N = 3 colours and Np massless quark flavours. 
2.2 NNLO infrared subtraction terms 

At NNLO, the m-jet production is induced by final states containing up to (to + 2) partons, 
including the one-loop virtual corrections to (to + l)-parton final states. As at NLO, one 
has to introduce subtraction terms for the (to + 1)- and (to + 2)-parton contributions. 
Schematically the NNLO m-jet cross section reads, 



VS,1 
'NNLO 



&<?NNLO — \ (do-NNLO ~ d(T NNLO 

) + / daf 
+ / i^NNLO ~ da NNLo) + / da l 

+ / d °NNLO > ( 2 - 14 ) 

where daff NLO denotes the real radiation subtraction term coinciding with the (m + 2)- 

parton tree level cross section dcr NNLO in all singular limits. Likewise, da^^ LO is the 

one- loop virtual subtraction term coinciding with the one-loop (m+ l)-parton cross section 

d a NN!LO m an singular limits. Finally, the two- loop correction to the m-parton cross section 
v 2 

is denoted by da^ NLQ . 
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2.3 Tree-level double real radiation subtraction terms 

Let us first consider the construction of the subtraction terms for the double radiation 
contribution dafj NLO , which shall correctly subtract all single and double unresolved sin- 
gularities contained in the (m + 2)-parton real radiation contribution to m-jet final states, 

d(T NNLO = N ^2 d $m+2(Pl, ■ ■ ■ ,Pm+2; Q)^— 
rn+2 

X \M m+2 { P U ■ ■ ■ ,Pm+2)\ 2 Jt +2 HP1, ■ ■ ■ ,Pm+2) • (2.15) 



S m +2 



Single real radiation singularities correspond to one parton becoming soft or collinear, while 
double real radiation singularities occur if two partons become soft or collinear simultane- 
ously. Singular terms in these limits can be identified by requiring a minimum number of 
invariants tending to zero in a given kinematical configuration. This number depends on 
the limit under consideration and follows from the phase space volume available to a given 
configuration. A detailed discussion of the kinematical definition of double unresolved 
limits can be found in [18, 19]. 

We must distinguish the following configurations: 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour connected but share a common radiator 
(almost colour- unconnected) . 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour- unconnected) . 

The first configuration was treated already in the context of antenna subtraction at 
NLO in Section 2.1 above. In the context of the construction of daf^ NLO , the same single- 
particle subtraction terms can be used. These do however not yet guarantee a finite 
(m + 2)-parton contribution in all single unresolved regions for two reasons: (1) while 
the jet function in da^ LO ensured that the subtraction term is non-zero only in the single 
unresolved limit it was constructed for, this is no longer the case for single unresolved 
radiation at NNLO; (2) the subtraction terms for the remaining three double unresolved 
configurations will in general be singular in the single unresolved regions, where they do 
not match the matrix element. Both problems will be addressed below. 

The remaining three configurations (b)-(d) are illustrated in Figures 4, 5 and 6. The 
singular behaviour of the full (m + 2)-parton matrix element in these configurations is the 
product of double unresolved factors (see Section 8 below) and reduced m-parton matrix 
elements. Subtraction terms for all these configurations can be constructed using either 
single four-parton antenna functions or products of two three-parton antenna functions. 
In all cases, attention has to be paid to the matching of different double and single unre- 
solved regions. This problem has been addressed already in several earlier publications on 
subtraction at NNLO [25-29,44], the most concise discussion can be found in [29]. 

In the following, we construct the subtraction terms for all four configurations. 
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Figure 4: Colour connection of the partons showing the parent and daughter partons for the 
double unresolved antenna. 
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Figure 5: Colour connection of the partons showing the parent and daughter partons for two 
adjacent single unresolved antennae. 
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Figure 6: Colour connection of the partons showing the parent and daughter partons for two 
disconnected single unresolved antennae. 



2.3.1 Subtraction terms for single unresolved partons 

Starting point for the subtraction terms for single unresolved partons are the NLO single 
unresolved antenna subtraction terms (2.5), 

1 



da 



S,a 

NNLO 



m+2 



d$ m+2 (pi, • • • ,p m +2;q)- 



~>m+2 



^2 X ijk \M m+ l(pi, ...,PI,PK,--- ,Pm+2)\ 2 4r +1) (Pl> • • • ,PI,PK,- ■ ■ ,Prn+2) 
j 

(2.16) 

where the NLO jet function is now replaced by J^ n+1 ' ) . In contrast to the NLO case, 
subtracting these terms from the full (m + 2)-parton matrix element does not ensure a finite 
contribution in all single unresolved regions. This behaviour is related to the fact that at 
NNLO the jet function J^ Ti+1 ' ) allows one of the (m + 1) momenta to become unresolved, 
while at NLO jffl required all m momenta to be hard. As a consequence, all momenta in 
the antenna function, including pj, could be resolved while one of the momenta present in 
the reduced (m + l)-parton matrix element becomes unresolved. We distinguish between 
two cases: (1) where pi or px become unresolved and (2) where any other momentum, 
p , in the matrix element becomes unresolved. Case (1) is necessarily a double unresolved 
limit since pi and px are linear combinations of two momenta; this configuration is treated 
below in the discussion of the double unresolved limits of (2.16). 

Case (2) is a single unresolved limit, since pj is resolved, while p a becomes unresolved. 



In this limit, da 



S,a 

NNLO 



becomes singular. Its singular structure in this limit does not coincide 



with the limit of the full (m + 2)-parton matrix element (which is already subtracted by 
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the appropriate antenna term in dcr^% LO when j = o). The unresolved momentum p Q in 
this situation is not colour-connected to the pj of X^- k . The form of the spurious singular 
terms in da^ a NLO is therefore equal to the product of the antenna function with a singular 
limit of the reduced (m + l)-parton matrix element. It is equivalent to the product of 
two almost colour-unconnected or colour-unconnected antenna functions with the reduced 
m-parton matrix element: 

X ijk \M m+ l(pi, . . . ,PI,PK, • • • ,Pn,Po,Pp, ■ ■ ■ ,Vm+2)\ 2 

„n „n , . . , _ _ _ _ \\2 P° unresolved 
- X ijk X nop \M m (pi, ■ ■ ■ ,PI,PK, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+2)\ ► 0, 

where Pn = Pk is allowed. Such structures appear as double unresolved subtraction terms 
for the configurations (c) and (d) defined above, and it will be shown in the appropriate 
subsections below, that the simple collinear limits of these double subtraction terms yield 
the correct behaviour to cancel the spurious single unresolved poles of du^ a NLO . 

To construct the double unresolved subtraction terms, we need to investigate the be- 
haviour of the subtraction term (2.16) in all its double unresolved limits. We distinguish 
the cases of colour-connected and colour-unconnected double unresolved limits. In the 
colour-connected, or type (b), double unresolved limits, (2.16) can be non- vanishing only 
for one particular momentum configuration: two neighbouring pairs of colour-connected 
momenta becoming independently collinear, where one of the pairs lies inside the antenna, 
while the other pair consists of the remaining antenna momentum and its colour-connected 
neighbour. Each configuration of this type is contained precisely twice in (2.16), since there 
are two possibilities of attributing the inside/outside pair. Where appropriate, we shall call 
these limits colour-neighbouring in the following, and decompose colour-connected limits 
into genuinely colour-connected and colour-neighbouring limits. In the colour-unconnected, 
type (c,d), double unresolved limits, (2.16) is always non- vanishing and yields twice the 
double unresolved limit of the (m + 2)-parton matrix element. This factor of two comes 
from the fact that each double unresolved limit requires both the unresolved antenna mo- 
mentum pj and one other momentum p Q appearing in the (m+ l)-parton matrix element to 
become unresolved; the role of the pj and p Q can be interchanged, resulting in two identical 
terms contributing to the same limit. 

To summarise the above discussion, we conclude that (2.16) yields twice the (m + 2)- 
parton matrix element in all colour-unconnected, almost colour-unconnected and colour- 
neighbouring double unresolved limits, while vanishing in all genuinely colour-connected 
double unresolved limits. 

2.3.2 Subtraction terms for two colour-connected unresolved partons 

When the two unresolved partons j and k are adjacent, we construct the subtraction 
term starting from the four-particle tree-level antenna X^ kl , which is an appropriately 
normalised four-particle matrix element [44,46,47]. By construction, it contains all colour- 
connected double unresolved limits of the (m + 2)-parton matrix element associated with 
partons j and k, but it can also be singular in single unresolved limits associated with 
j or k, where it does not coincide with limits of the matrix element. To ensure a finite 
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Figure 7: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles are 
colour connected. The term in square brackets represents both the antenna function Xfj kl and the 
antenna phase space d&x ijkl • 



subtraction term in all these single unresolved limits, we therefore subtract the appropriate 
limits of the four-particle tree antennae, which are products of two tree-level three-particle 
antennae, such that the colour-connected double subtraction term reads: 



dcr NNLO = N ^2 d$m+2(pi, . . . ,p m +2] Q)-q— 



m+2 Jm + 2 



[.^■ijkl - ^ijk^IKl - ^jkl^iJL) 

jk 



\M m (pi, . . . ,PI,PL, ■ ■ ■ ,Pm+2)\ 2 J^(Pl, • • • ,PI,PL, • • • ,Pm+2) 



(2.17) 



where the sum runs over all colour- adjacent pairs j, k and implies the appropriate selection 
of hard momenta 

As before, the subtraction term involves the m-parton amplitude evaluated with on- 
shell momenta pi, ■ ■ ■ ,Pi,Pl, ■ ■ ■ ,Pm+2 where now pi and pi are a linear combination of pi, 
Pj, pk and p\. As for the NLO antenna of the previous section, the tree antenna function 
X®j kl depends only on Pi,Pj,Pk,Pi- Particles i and I play the role of the radiators while 
j and k are the radiated partons. Antenna factorisation of squared matrix element and 
phase space in this configuration is illustrated pictorially in Figure 7. 

Once again, the jet function in (2.17) depends only on the parent momenta p~i,p~L 
and not Pi, ■ ■ ■ ,pi- One can therefore carry out the integration over the unresolved antenna 
phase space (or part thereof) analytically, exploiting the factorisation of the phase space, 

d$ m+2 (pi,...,p m+ 2;g) = d^mip!, . . . ,pj,p L , . . . ,p m+1 ;q) ■ d<S>x ijkl (Pi,Pj,Pk,Pf,Pi + Pl) ■ 

(2.18) 

The factorisation [21,25,31] is obtained by redefining a set of four massless on-shell momenta 
(radiator, two unresolved partons, radiator) into two on-shell massless momenta. The 
NNLO antenna phase space d$x ijM is proportional to the four-particle phase space, as can 
be seen by using m = 2 in the above formula such that 

d*4 = P 2 d* Xijkl . (2.19) 
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One suitable representation of momenta for analytic integration purposes is the tripole 
phase space mapping of [31] where parton i is the emitter and parton I the spectator, 

f'f^^-* «-t4s* (2 ' 20) 

with 

Pi-Pi + Pl-Pj + Pl-Pk , 00 ,x 

= : ; ; : ; • K 2 - 21 ) 

Pi-Pj+ pi .p k + Pi -Pi + Pj -Pk + Pj -Pi + Pk -Pi 

For numerical implementation the mapping of [25] is suitable. 

For the analytic integration, we can use (2.18) to rewrite each of the genuine four- 
particle subtraction terms in the form, 

\M m \ 2 Jt ] d$ m J &$x ljkl X% kh (2.22) 

as before. The analytic integral of the four-particle antenna subtraction term is there- 
fore the antenna function integrated over the fully inclusive antenna phase space, again 
including a normalisation factor to account for powers of the QCD coupling constant, 

= ^ e en ) 2 j d ^ x m ■ (2-23) 

This integration is performed analytically in d dimensions using the techniques explained 
in [31] to make the infrared singularities explicit. 

The double unresolved subtraction term (2.17) correctly approximates the (m + 2)- 
parton matrix element (2.15) in all double unresolved limits where colour-connected partons 
j and k become unresolved, except in the configuration where parton j becomes collinear 
with parton i, while parton k becomes collinear with parton I. In this (so-called colour- 
neighbouring) configuration, the tree-level four-parton antenna function Xijki correctly 
approximates the (m + 2)-parton matrix element, while the products of tree-level three- 
parton antenna functions X^ k Xj Kl and Xj kl Xf JL in (2.17) are non- vanishing, and are each 
equal to the double unresolved limit of the matrix element. This colour-neighbouring mo- 
mentum configuration was already discussed in the context of the double unresolved limits 
of d<r}^ LO in the previous subsection. It can be seen easily that the two products of tree- 
level three-parton antenna functions X^ k Xj Kl and Xj M Xf JL in (2.17) exactly cancel the 
spurious double unresolved singular terms in present in (2.16) in these colour-neighbouring 
double unresolved limits. 

In all genuinely colour-connected limits, the four-parton tree-level antenna function 
X^ki in (2.17) correctly matches the singular structure of the (m+2)-parton matrix element, 
while the products of tree-level three-parton antenna functions Xfj k Xj Kl and X® kl X® JL 
vanish (although these terms are singular at first sight, they are not sufficiently singular to 
overcome the small phase space volume associated with the double unresolved limits under 
consideration). 

By construction, (2.17) vanishes in all single unresolved limits: the antenna functions 
cancel each other in these limits, and singularities in the m-parton matrix element are 
forbidden by the jet function. 
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Figure 8: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m' + 2)-particle phase space when the unresolved particles are 
almost colour-unconnected. 



2.3.3 Subtraction terms for two almost colour-unconnected unresolved partons 

There are double unresolved configurations where the unresolved partons are separated by 
a hard radiator parton, for example, i,j,k,l,m where j and I are unresolved as illustrated 
in Fig. 8. These configurations are called almost colour- unconnected. In this case, we take 
the strongly ordered approach [41] where i,j, k form an antenna with hard partons I and 
K yielding an ordered amplitude involving I,K,l,m. As usual, the momenta of the hard 
radiator partons I and K are labelled pi and pk and are constructed from Pi,Pj,Pk- The 
cases where I is unresolved are then treated using an antenna K,l,m with hard partons 
K and L. Here the momenta of the hard radiator partons K and L are labelled px and 
PM which are made from pK,Pi,Pm- The other case where first k,l,m form an antenna 
followed by i,j, K is also included. 

In constructing the subtraction term, we have to account for the fact that the previ- 
ously defined subtraction term d<7^ i0 (2.16) already contributes in the almost colour- 
unconnected configurations by subtracting twice the limit of the {m! + 2)-parton matrix 
element in these limits (to avoid confusion between the indivial unresolved momenta and 
the total number of particles, we denote the latter {m! + 2) instead of (m + 2) in this 
subsection). Therefore, the almost colour- unconnected subtraction term reads, 

d(J NNLO = d< ^m'+2(Pl, ■ ■ ■ ,Pm'+2;q)-^— 

m'+2 Sm '+ 2 

X ijk X ™lK \ M m'{Pl, ■ ■ ■ ,PI,PK,PM, ■ ■ ■ ,Prn'+2)\ 2 

hi 

XJ^ \P1, • • • ,PI,PK,PM, ■ ■ ■ ,Pm'+2) 

+ X klm x ijK \M m '(pi, ■ ■ ■ ,Pl,PK,PM, • • • ,Pm'+2)\ 2 

hi 

T (m') I » »s \ 

XJrn* (Pl> • • • &I>PK,PM, ■ ■ ■ ,Pm'+2) 



(2.24) 



where x® mlK denotes a sub-antenna, containing only the collinear limit of m with /, but 
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Figure 9: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles are 
colour- unconnected. 



not the collinear limit of I with K; in the soft limit of I, this sub-antenna yields half the 

(m') 

soft eikonal factor. The jet function J; m , and the matrix element depend only on the 
momenta associated with the m'-particle final state. Therefore, the integral over the phase 
space associated with the leftmost antenna function can be carried out analytically, like in 
the single unresolved case discussed above. After this integration, the sub-antennae can be 
recombined to full antenna functions again. 

If required, it is also possible to carry out the integrals over both phase spaces associ- 
ated with the two antenna functions. In this case, one can exploit the factorisation of the 
phase space to analytically perform the integrals, 

d$m'+ 2 (Pl, • • • ,Pm'+2]q) = 

d$m'+l(Pl, • • • ,PI,PK, • • • ,Pm'+2;q) " ^X ljk (Pi,Pj , Pk) 
= d$ m /(pi, ...,PI,PK,PM,--- ,Pm'+2; q) ■ d<f>x Klm (PK,Pl,P m ;PK +Pm) 

■d®x ijk (Pi,Pj,Pk;Pi +Pk) ■ (2.25) 

In carrying out the analytic integration, some care has to be taken in the second dipole in- 
tegral d&x Klm , which will pick up e-dependent factors from the first integral (both integrals 
are fully independent only in four dimensions) related to the correlation of unit volumes 
involved in the two integrals. As a consequence, the analytic integration of the above 
formula will not yield the product of two independent integrated NLO antenna functions. 

Concerning the single unresolved limits of da^f NLO , (2.24) becomes singular only if 
one of the antenna partons j or I is unresolved. In these cases, it exactly cancels the 
spurious single unresolved singularities encountered in da S ^ a NLO for the configuration of an 
unresolved momentum p Q in the {m! + l)-parton matrix element, which is almost colour- 
unconnected to the momenta in the antenna function. In the double unresolved limits of 
d<rfjx L Q, momenta j and / are both unresolved, yielding minus the double unresolved limit 
of the (m' + 2)-parton matrix element, as required to cancel the oversubtraction incurred 
by da^ a NLO in these limits. 

2.3.4 Subtraction terms for two colour-unconnected unresolved partons 

When two unresolved partons j and o are completely disconnected, the (m + 2)-parton 
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matrix element factorises into the product of two uncorrelated single unresolved factors 
with a reduced m-parton matrix element. In this limit, the single unresolved subtraction 
term da^ LO (2.16) subtracts a contribution which is twice the limit of the full (m + 2)- 
parton matrix element. The dedicated colour- unconnected double unresolved subtraction 
term has to compensate for this oversubtraction. This subtraction term is a sum over 
independent three-particle tree-level antennae as illustrated in Figure 9, 



d(7 NNLO = ^2 d< ^rn+2(pi,- ■ ■ ,p m +2;q)-^ 



m+2 



m+2 



^2 X ijk X nop \M m (pi, • • • ,PI,PK, • • • ,PN,PP, • • • ,Pm+2) 



3,0 



J m\pii ■ ■ ■ ,PI,PK, • • • ,PN,PP, • • • ,Pm+2) 



(2.26) 



where the summation over o is such that it only includes colour-connected antenna config- 
urations X® op which have no common momentum with X^- k . 

The subtraction term involves the m-parton amplitude evaluated with on-shell mo- 
menta pi, ■ ■ ■ ,Pi,Pk, ■ ■ ■ 1PN1PP1 ■ ■ ■ ,Pm+2 where now pj and px are a linear combination 
of pi, pj, pk and pN and pp are made from p n , p Q and p p . The sum in the above equation 
is such that no product of two antenna configurations appears twice. Each antenna is fully 
independent, and the phase space factorises as, 

d$ m+2 (pi, • • • ,Prn+2; q) = d$ m (pi, . . . ,Pl,pK, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+2] q) 

x ^x lJk {puPj,Pk;Pi + Pk) ■ d$x nop (Pn,Po,p P ;PN +pp) ■ 

(2.27) 

This contribution is either integrated only once, and then cancelled with subtraction terms 
present in the (m + l)-parton contribution, or integrated twice to make its full infrared 
pole structure explicit. Integration can be obtained by repeated use of (2.11), so that each 
term in the sum is in the form, 

|-M m | J^j ^ijki s ijk) X nop (s n0 p). 

The subtraction term (2.26) contributes to single unresolved limits if either pj or p D 
are unresolved. In these limits, it cancels precisely the spurious poles associated with the 
momentum p a becoming unresolved in the (m + l)-parton matrix element of (2.16). In 
the double unresolved limits (pj and p a ) it is designed for, it cancels the oversubtraction 
incurred by (2.16), thus ensuring a proper subtraction of the singular terms of the (m + 2)- 
parton matrix element in these. 

2.3.5 Correction terms in the m-jet region 

The full double radiation subtraction term is given as sum of all subtraction terms con- 
structed above: 

d&NNLO = da NNLO + da NNLO + da NNLO + da NNLO ■ (2.28) 
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As outlined in the previous subsections, this subtraction term correctly approximates the 
(m + 2)-parton matrix element contribution to m-jet final states as defined in (2.15) in all 
double and single unresolved regions. Although individual terms in (2.28) contain spurious 
singularities in these limits, they cancel among each other in the sum. 

The integrated form of (a) corresponds to an (m + l)-parton configuration, while the 
integrated forms of (b), (c) and (d) are either (m + l)-parton or m-parton configurations 
(for all but the four-parton antenna terms in (b), we can actually choose which type of 
configuration we want to integrate). They are added with the two- loop m-parton and 
the one- loop (m + l)-parton contributions to m-jet final states to yield integrands free of 
explicit infrared poles. 

2.4 Single unresolved loop subtraction terms 

The (m + l)-parton one-loop contribution to m-jet final states at NNLO is 



where we introduced a shorthand notation for the one-loop corrected contribution to the 
(m + l)-parton squared matrix element, 



• • • ,Pm+i)\ 2 = 2Re (AC+ibi, • • • ,Pm+i)*0*(pi, . . . , Pm+ i)) , (2.30) 



where the sub leading contributions in colour are again implicit. This expression contains 
two types of infrared singularities. The renormalised one-loop virtual correction to 
the (m + l)-parton matrix contains explicit infrared poles, which can be expressed using 
the infrared singularity operators defined in Section 4 below. On the other hand, one of the 
(m + l)-partons can be unresolved, leading to infrared singularities which become explicit 
only after integration over the unresolved patch of the final state (m + l)-parton phase 
space. The single unresolved limits of one- and two-loop amplitudes are investigated in 
detail in [11-15]. 

To carry out the numerical integration over the (m + l)-parton phase, weighted by 

VS 1 

the appropriate jet function, we have to construct an infrared subtraction term da NI J LO , 
which purpose is twofold: on the one hand, it should fully account for the explicit infrared 
singularities arising from the one-loop correction, while on the other hand it should approx- 
imate (2.29) in all single unresolved limits. After subtraction of <ia^f^ L Q from da^l^ LO , 
the integrand itself should be free from explicit infrared poles, and all singular limits due to 
single unresolved phase space configurations must be properly accounted for. We require 
three types of subtraction terms: 

(a) The explicit infrared poles of the (m+l)-parton one-loop matrix element are removed. 

(b) The single unresolved limits of the (m + l)-parton one-loop matrix element are sub- 



da 



.v,i 

NNLO 




x|^m+l(Pl,--->Pm+l)| 2 4T +1) (P1,...,IWi) , 



(2.29) 



tracted. 



(c) Oversubtracted explicit and implicit pole terms are removed. 
All three types of subtraction terms are constructed below. 
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2.4.1 Subtraction of explicit infrared poles 

It is a well known fact from NLO calculations, that the explicit infrared poles of one-loop 
matrix elements cancel with the corresponding infrared poles obtained by integrating out 
all single unresolved configurations from the real radiation matrix elements contributing 
to the same (infrared safe) observable [49]. We can therefore subtract all explicit poles 
present in da^^ LO with the subtraction term 



da NNLo = N d ^ >m + 1 (pi ' • • • > Pm+i ; q) ~q~ 



Sm+1 



ik 



~ X ijk{ s ik) \M m +l(Pl, • • • ,Pi,Pk, • • • ,Pm+l)\ 2 Jm +1) iPli ■ ■ ■ ,Pi,Pk, • • • ,Pm+l) 

(2.31 



where the sum runs over all colour-connected pairs of momenta (pi,Pk)- The symbol 
^■ijk( s ik) denotes the normalised integral of the three-parton antenna function, with a 
parton of type j emitted between partons i and k, with the total invariant mass of the 
antenna particles given by the invariant mass of pi and pk- Since in the single unresolved 
mapping of the antenna phase space (2.8), pi + pj + pu = Pi +Pk, we can identify the 
momenta Pi and pk of the (m+ l)-parton phase space with the momenta fij and px obtained 
in the factorisation of the (m + 2)-parton phase space in Section 2.3.1. Consequently, using 
the antenna factorisation of the (m + 2)-parton phase space in (2.16) and carrying out the 
phase space integral for each antenna X^, we find 

da NNLO = ~ da NNLO ' (2.32) 

which is of course merely a consequence of the cancellation of infrared poles at NLO. 

It is evident that subtraction of dcrJ^^Q from dcr]^ i0 ensures an integrand free 
from explicit infrared poles over the whole region of integration. In the single unresolved 
regions, both these functions develop further infrared singularities, which do not coincide. 
Therefore, we have to introduce further subtraction terms to ensure a finite integrand in 
all single unresolved regions. 

2.4.2 Subtraction terms for one-loop single-unresolved contributions 

In simple unresolved limits, the behaviour of (m + l)-parton one-loop amplitudes is de- 
scribed by the sum of two different contributions [11-15]: a simple unresolved tree level 
factor times a m-parton one-loop amplitude and a simple unresolved one-loop factor times 
a m-parton tree-level amplitude, as illustrated in Figure 10. Accordingly, we construct the 
one-loop single unresolved subtraction term as 

da NNLO = ^ Yl d $m+l(Pl, ■ ■ ■ ,p m+ i;q)- 



m+1 Jm+1 



E 



X iik iMUPli • • • ,PliPK, ■ ■ ■ ,Pm+l)| 2 4r } (Pl> • • • ,PI,PK, • • • ,Pm+l) 
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Figure 10: Illustration of NNLO antenna factorisation representing the factorisation of both the 
one- loop "squared" matrix elements (represented by the white blob) and the (m+ l)-particle phase 
space when the unresolved particles are colour connected. The terms in square brackets represent 
both the three-particle tree-level antenna function X^- k and the three-particle one-loop antenna 
function Xj- k and the antenna phase space. 

+ X ijk \M m (Pl, ...,PI,PK,--- ,Pm+l)\ 2 ^{Pl, ■ ■ ■ ,Pl,PK, ■ ■ ■ ,Pm+l) ■ 

(2.33) 

In this expression, we have introduced the one-loop three-parton antenna function X^ k , 
which depends only on the antenna momenta Pi,Pj,Pk- It correctly describes all simple 
unresolved limits of the difference between an (m + l)-parton one- loop corrected squared 
matrix element and the product of a tree-level antenna function with the m-parton one- 
loop corrected squared matrix element. It can therefore be constructed out of one-loop 
three-parton and two-parton matrix elements, as outlined in Section 3 below. It should 
be noted that X\- k is renormalised at a scale corresponding to the invariant mass of the 
antenna partons, Sij k , while the one-loop (m + l)-parton matrix element is renormalised at 
a scale fi 2 (which is often chosen to be q 2 in jet production in e + e~-collisions). To ensure 
correct subtraction of terms arising from renormalisation, we have to substitute 

Xl jk - Xl jk + ^ ^pAg- fc ((s ljk r - (» 2 r) (2-34) 

in (2.33). The terms arising from this substitution will in general be kept apart in the 
construction of the colour-ordered subtraction terms, since they all share a common colour 
structure (3q- 

In contrast to all other expressions appearing in the construction of the one-loop single 
unresolved subtraction terms, Xj- k can never be related to integrals of tree-level subtraction 
terms. Therefore, this component of the subtraction term must cancel with parts of the 
two-loop m-parton contribution and we must integrate it over the three-parton antenna 
phase space. This can be accomplished using the techniques described in [31] and yields 

4*( a 0"*) = (Sir 2 (4vrr ^) j d$ Xijk X} jk . (2.35) 
2.4.3 Compensation terms for oversubtracted poles 

By construction, (2.33) correctly approximates the one-loop (m + l)-parton contribution to 
m-jet final states in all single unresolved limits. However, outside these limits, where (m + 
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l)-parton configurations form m-jet final states because of the experimental jet resolution 
criteria, this subtraction term no longer coincides with the squared matrix element, and 
induces spurious explicit infrared poles. 

To compensate for these explicit infrared poles outside the singular regions, we intro- 
duce a further subtraction term, 

d(T NNLO = N d ®rn+l (pi , • • • , Pm+1 5 q) 

,Pi,Pk, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+l)\ 



ik 



x4r ) (pi.--->p»> 



Pfc, - ■• ,PN,PP, ■ ■ ■ ,Pm+lj 



(2.36) 



where any of the two momenta pi and pk can be equal to p n , p p , pn or ftp, but not to 
the unresolved momentum p a . Individual terms in do^fj^Q cancel with individual terms 
in da^ b NLO (2.17), daf} c NLO (2.24) and da S jf NLO (2.26), after one of the two three-parton 
antenna phase space integrals is carried out. Any remaining terms of da^^Q can be 
integrated over the antenna phase space d&x nop ■ In the case of pi or p^ coinciding with 
one of the momenta of the antenna phase space, some care has to be taken in carrying out 
the integrals, which differ from the standard tree-level three-parton antenna integrals by 
normalisation factors coming from X^ k (sik)- 

The subtraction term (2.36) cancels the explicit infrared poles of (2.33) in the region 
where all (m+1) partons are theoretically resolved, thus ensuring a finite integrand. In the 
single unresolved regions, the sum of (2.36) and (2.31) vanishes, as can be seen rather easily 
from the nature of NLO antenna subtraction. Consequently, only da^ S ^ contributes in 
this region, as required to cancel the singularities of the (m + l)-parton squared matrix 
element. 

2.4.4 Correction terms in the m-jet region 

The full one-loop single real radiation subtraction term is the sum of all subtraction terms 
constructed above: 

i VS,1 _ i VS,l,a , i VS,l,b , VS,l,c (0 o 7 \ 

aa NNLO ~ aa NNLO aa NNLO aa NNLO m \ z - 6 <) 

As outlined in the previous subsections, this subtraction term correctly approximates the 
one-loop (m+l)-parton squared matrix element contribution to m-jet final states as defined 
in (2.29) in all single unresolved regions and removes all explicit infrared poles. 

While da^f^^Q cancels fully with da^ a NLO , parts of the remaining two terms have to 
be integrated to yield m-parton configurations, which are then added with the two-loop 
m-parton contributions da^ NLQ . The precise nature of cancellations between the terms 
appearing in the integrated forms of da^pjLo anc ^ terms in the unintegrated da^^Q^ 
differs considerably among the different colour structures. Therefore, no generic formula 
for these cancellations can be stated. In Section 9, we will illustrate these cancellations on 
the example of the subleading colour contribution to e + e~ — > 3j at NNLO. 
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2.5 Comparison with other approaches 

The NNLO antenna subtraction approach which we derive in detail in this section was previ- 
ously sketched in our earlier publications [44,45], where specific applications to e + e~ — > 2j 
and e + e~ — > 3j were considered. The case two-jet production is special in several as- 
pects, particularly since the double unresolved tree-level and single unresolved one-loop 
subtraction terms exactly coincide with the full matrix elements they are supposed to ap- 
proximate. Moreover, the jet functions for two- and three-parton final states fulfil special 
relations which are in general not fulfilled for m- and (m + l)-parton final states. Most 
notably (see (5) in [45]), 

(3) (2) (3) 

J 3 (PhPj,Pk) = J-2 (Pi,Pk) ~ J<2 {Pi,Pj,Pk) , (2.38) 

which can be used to express (3.17) of [44] in the form of (2.28) presented in this paper. 

Several other approaches to the handling of infrared singularities at NNLO have been 
proposed in the literature. One can distinguish two substantially different lines of making 
the infrared singularities in real radiation at NNLO explicit: either through a direct expan- 
sion of matrix element and phase space or through subtraction terms. In expanding the 
full real radiation matrix elements and the full multi-parton phase space in dimensional 
regularisation through sector decomposition [30,31], one arrives at a Laurent expansion in 
e with coefficients in terms of distributions in the Lorentz invariants associated with the 
process. These coefficients still contain the full phase space integration, but they are finite 
and the integrals can be carried out numerically. Using this approach, NNLO results have 
been obtained for e + e~ — ► 2j [32], pp —> H + X [33] and more recently muon decay [34] . In 
approaches invoking subtraction terms, such as the antenna subtraction proposed here, the 
real radiation singularities from multiparton matrix elements are removed by subtracting 
appropriate approximations, which coincide with the matrix elements in all singular limits, 
but are sufficiently simple to be integrated analytically. 

A first formulation of the antenna subtraction method at NNLO was presented in 
Refs. [13,25]. In this work, all antenna phase space mappings required at NNLO were 
derived and documented. These mappings are to be used in the implementation of our 
formulation of NNLO antenna subtraction as well. Concerning the construction of sub- 
traction terms, Ref. [25] addresses processes involving gluons only, and restricts its results 
to the leading colour contributions (which are sufficient in purely gluonic processes at 
NNLO). The antenna subtraction terms in [25] are presented only in their unintegrated 
form. To implement them in an actual calculation, they still need to be integrated over 
the appropriate antenna phase spaces to make their infrared singularities explicit. 

A generalisation of the NLO dipole subtraction formalism to NNLO was presented 
in [14,26], where subtraction terms for all partonic configurations in all colour structures 
were derived. Again, these were presented only in their unintegrated form, while their 
analytical integration over the NNLO dipole phase spaces is a necessary prerequisite for 
their implementation into a numerical programme. 

Finally, a third approach to the construction of subtraction terms at NNLO was pre- 
sented in [28] , where an iterative procedure of subtraction and subsequent cancellation of 



- 22 - 



oversubtracted terms is invoked. Using this procedure, NNLO subtraction terms for the 
CfTr colour factor in e + e~ — ► 2j are constructed in [28], where they are also integrated 
analytically and implemented into a numerical programme. The agreement of [28] with 
earlier results [32,44] illustrates the potential of this method. However, the full set of sub- 
traction terms for all relevant partonic configurations in all colour structures is at present 
not yet derived (and thus not available in integrated form either) for this method. 

More specific issues about subtraction at NNLO were addressed in two further works: 
[27] deals with overlapping collinear divergences in the initial state, and proposes a method 
for their systematic separation. More recently, a detailed study of the matching of single 
and double unresolved regions in subtraction at NNLO was presented in [29]. 

To implement a subtraction method at NNLO, one requires explicit expressions for all 
subtraction terms (double real radiation at tree level and single radiation at one-loop) in 
their unintegrated form, as well as the integrals of these expressions over the appropriate 
(double unresolved or single unresolved) phase spaces. Up to now, none of the approaches 
available in the literature provided these to an extent which would permit implementation 
for a realistic observable. 

In the following sections, we will present the NLO and NNLO antenna functions for all 
partonic configurations and all colour factors. These are the building blocks of NLO and 
NNLO antenna subtraction terms required to implement our formulation of the colour- 
ordered antenna subtraction method. We provide the integrals of these subtraction terms 
for the kinematical situation of final state radiation. The method itself can in principle 
be extended to processes with hadrons in the initial state. In this case, the antenna 
subtraction terms could be constructed from the same building blocks, but different phase 
space factorisations and consequently different integrals are required. 

3. Notation and structure of antenna functions 

In this paper we derive antennae for all possible pairs of hard partons, quark-antiquark 
(Section 5), quark-gluon (Section 6) and gluon-gluon (Section 7). Underlying the antenna 
is a colour connected pair of hard partons, that emit radiation between them. The three- 
particle antennae involve one unresolved parton, while the four-particle antennae involve 
two unresolved partons. The antenna may be at the tree or one-loop level. 

Each antenna is determined by both the external state and the pair of hard partons it 
collapses to. In general we denote the antenna function as X. For antennae that collapse 
onto a hard quark-antiquark pair, X = A for qgq and qggq, X = B for qq'q'q and X = C for 
qqqq final states. Similarly, for quark-gluon antenna, we have X = D for qgg and qggg and 
X = E for qq'q' and qq'q'g final states. Finally, we characterise the gluon-gluon antennae 
as X = F for ggg and gggg, X = G for gqq and gqqg and X = H for qqq'q' final states. 
Some of these antenna functions decompose into a leading colour and a subleading colour 
contribution. Where appropriate, the subleading colour contribution to X is denoted as 
X. Finally, the one loop antenna functions also contain contributions from closed quark 
loops, which we denote as X . The notation for the different tree-level and one- loop antenna 
functions is summarised in Table 1. 
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tree level 


one loop 


quark-antiquark 






qgq 


A° 3 {q,g,q) 


Al(q,g,q), Al(q,g,q), A\{q,g,q) 


qggq 


Al(q,g,g,q), Al(q,g,g,q) 




qq'q'q 


B%(q,q',q',q) 




qqqq 


C2(q,q,q, q) 




quark-gluon 






qgg 


D° 3 (q,g ig ) 


Dl(q,g,g), D\{q,g,g) 


qggg 


Dl(q,g,g,g) 




qq'q' 


El(q,q',q') 


El(q,q',q'), E\{q,4,tf), El(q,q',q') 


qq'q'g 


El{q,q',q',g), E° 4 (q,q' ,q> , g) 




gluon-gluon 






ggg 


F$(g,g,g) 


F 3 {g,g,g), F^{g,g,g) 


gggg 


F%(g,g,g,g) 




gqq 


G° 3 (g,q,q) 


Gl(g,q,q), Gl(g,q,q), G\{g,q,q) 


gqqg 


Gl(g,q,q,g), Gl(g,q,q,g) 




qqq'q' 


H2(q,q,q',q') 





Table 1: List of tree- level and one-loop colour-ordered antenna functions. The tilde denotes 
subleading colour contributions, and the hat flavour-number dependent corrections. 

Each antenna is a function of the invariants formed by the momenta of the final state 
particles. 

It will prove useful to introduce the operators Voles and J-inite that select either the 
singular or finite contribution from a particular antenna. For example, 

Voles(X) 

extracts the singular contribution from the antenna X after integration over the unresolved 
phase space in terms of the infrared singularity operators of Section 4. It should be noted 
that Voles(X) also contains finite terms arising from the expansion of the infrared singular- 
ity operators. These terms are on the one hand of the type (elns) n (where s is the invariant 
mass of any pair of particles contained in the antenna, or the total invariant mass of all 
antenna particles) , on the other hand they contain transcendental constants resulting from 
the expansion of the normalisation factors appearing in the different infrared singularity 
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operators. To extract the remaining finite contribution, we introduce 



Tinite(X) = X - Voles(X) . 

Generally, 

X = Voles(X) + Finite(X) + 0(e). (3.1) 

The one-loop antenna functions contain explicit poles from the loop integration. There- 
fore, the operators Voles and Finite can also be applied to their unintegrated forms X. 
The action of these operators is again to decompose the unintegrated antenna in terms of 
infrared singularity operators describing the pole terms and a finite remainder. 

All antenna functions are derived from physical matrix elements: the quark-antiquark 
antenna functions from 7* — ► qq + (partons) [44] , the quark-gluon antenna functions from 
X —>■ g + (partons) [46] and the gluon-gluon antenna functions from H — ► (partons) [47]. 
The tree-level antenna functions are obtained by normalising the colour-ordered three- and 
four-parton tree-level squared matrix elements to the squared matrix element for the basic 
two-parton process, 

^ijk — D ijk,IK |^0 |2 ' 



I-Kli 



A iji _ b ijkl,IL TTltTTT ' \ 6 - 1 ) 



where S denotes the symmetry factor associated to the antenna, which accounts both for 
potential identical particle symmetries and for the presence of more than one antenna 
in the basic two-parton process. The one-loop antenna functions are obtained from the 
colour-ordered renormalised one-loop three-parton matrix elements as 

\M l -A 2 lA/f 1 i 2 

yl o 1 l^fcj Y \ JVI IK\ (o o\ 

X ijk - b ijkJK - X ijk j^-^ . [6.6) 

The numerical implementation of the three- and four-parton antenna phase space [25] 
requires the partonic emissions to be ordered. Ordering of emissions means that the two 
hard radiator partons defining the antenna are identified, and that each unresolved par- 
ton can become singular only with the two particles which are adjacent to it, i.e. with 
the two radiators for three-parton antenna functions and with one radiator and with the 
other unresolved parton for the four-parton antenna functions. For the sake of numeri- 
cal implementation, this implies two requirements: (1) the separation of multiple antenna 
configurations present in a single antenna function for three- and four-parton antenna func- 
tions and (2) the separation of non-ordered emissions (present only at subleading colour 
in the four-parton antenna functions) into terms that can be identified with a particular 
ordering of the momenta. 

In the colour-ordered quark-gluon and gluon-gluon antenna functions derived from 
physical matrix elements for neutralino decay [46] and Higgs boson decay [47], it is in 
general not possible to identify the hard radiators and the unresolved partons in a unique 
manner. The reason for this ambiguity is in the cyclic nature of the colour orderings, which 
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becomes evident already in the three-parton antenna functions: each pair of two partons 
can in principle act as hard radiators, resulting in more than one antenna configuration 
present in a single antenna function. For the three-parton antenna functions D® and F®, 
we illustrate in (6.13) and (7.13) below how these different antenna configurations can 
be disentangled, resulting in new sub-antenna functions where the hard radiators can be 
uniquely identified. Such a decomposition is also possible for four-parton antenna functions 
which display the same ambiguity. However, in the four-parton case repeated partial frac- 
tioning is required to extract individual sub-antenna configurations. Since this procedure 
introduces new denominators in the antenna terms, we will use it only for the numerical im- 
plementation. At this point, it should be pointed out that the cyclic ambiguity is inevitable 
if physical matrix elements involving gluons as radiators are used to derive the antenna 
functions, since the cyclic symmetry in colour space is enforced by gauge invariance. 

The decomposition of non-ordered emissions into different terms is discussed in detail 
in Section 5 using the quark-antiquark antenna function A®, which describes the emission 
of two gluons at subleading colour, as an example. In this antenna the gluons behave 
effectively like photons, coupling only to the quark and the antiquark. Consequently, each 
gluon can become collinear with both hard radiators. In (5.28) below, we illustrate how this 
antenna function can be decomposed into individual terms corresponding to different well- 
defined orderings by repeated partial fractioning in the invariants. The same procedure has 
to be applied for all other antenna functions which are not sufficiently ordered. Since this 
repeated partial fractioning yields increases the number of terms in the analytic expressions 
for the antenna functions quite considerably, we restrict ourselves to quoting only the full 
antenna functions, not their ordered decompositions. 

Repeated partial fractioning generates denominators which do not correspond to phys- 
ical propagators of four-parton matrix elements. Therefore, it is in general not possible 
to integrate the ordered terms analytically, at least not by merely applying the methods 
of [31]. Analytical phase space integration of functions involving non-propagator-type de- 
nominators is in principle possible [15,16] and requires a more involved reduction procedure 
and a larger set of master integrals. To avoid this problem we use the ordered sub-antenna 
functions in the numerical implementation, and ensure that all ordered contributions to 
a given antenna function are taken together with the same phase space factorisation (but 
different phase space mappings). With this procedure, the ordered contributions to the 
sub-antenna functions can be recombined to form the full antenna functions (related to 
physical matrix elements) quoted in the following three sections, which are then integrated 
analytically to make their infrared pole structure explicit. This recombination may be 
obtained by considering more than one colour-ordering for the full (m + 2)-parton matrix 
element, and by explicit subtraction of unwanted spurious singularities. This problem is 
discussed more explicitly using the four-parton quark-gluon antenna function D® of Sec- 
tions 6 and 8.3.2 as an example. 

4. Colour-ordered infrared singularity operators 

In order to express the singularity structure in a way that cancellations between real ra- 
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diation and virtual radiation can be made explicit, it is convenient to extract the infrared 
singularity structure using the I*- 1 ) -operator [48]. This operator describes the singular- 
ity structure of virtual one-loop amplitudes. At the two-loop level, the infrared singularity 
structure is described by -operators, a hard radiation operator H^ 2 ), the QCD /3-function 
and the gluon splitting constant K. The same pole structure is recovered [44,46,47] in the 
sum of double real radiation and one-loop single real radiation contributions, such that the 
sum of double-virtual, virtual single-unresolved and double-unresolved corrections is finite. 
Some cancellations, related to the one-loop correction to the soft gluon current S^ 2 ), take 
place only between virtual single-unresolved and double-unresolved corrections. 

In the formulation of [48], iW is a tensor in colour space, and contains imaginary parts 
from the analytic continuation of loop amplitudes from the Euclidian to the Minkowskian 
region. In the following, we will only consider colour-ordered matrix elements, for which 
iW is a scalar in colour space. Moreover, using the I^-operator to describe real radiation 
singularities, only its real part is relevant, since contributions from its imaginary part 
cancel once the double-virtual two-loop-times-tree-level and one-loop-times-one-loop parts 
are added together. 

The real radiation infrared singularity operators that appear in the integrated form of 
an antenna function are denoted by (ij = qq,qg, gq, gg) 



r(i) 



and 



(e s-) 



where the second operator describes contributions arising from the splitting of a gluon into 
a quark-antiquark pair and which are proportional to the number of light quark flavours 
N F . 



The operators are given by 

L qq l fc ) b qq) " 



^1 



^-qg( e i s qg) 



I (1) fe s ) 

gg v ' gg> 



2r(l - e) 
2r(l-e) 
'2r(l-e) 



e 2 + 2e 



e 2 + 3e 
1 11 



Re(-s q q)~ 



Re(— s, 



'in i 



6e 



Re(-s 



.'/// 1 



(e, S q g) = 



^qg,F( e i s qg) 



=«7 



2T 1 



l gg,F 



( e > s gg) 



=«7 



1 

6e 
1 



2T(1 - e) 3e 



Re(-s qg y 
Re(-s gg Y 



(4.1) 



The antiquark-gluon operators are obtained by charge conjugation: 



I (1) (e s 

*-qq l fc ) 6 



A <7.9 v 



e, s 



and 



r(l) 



.(e,s 



g<i i 



^qg,F( e i s gq) 



"■gq v-i a gqi ~ *-qg \^^ J gq> gq,F\ 

At NNLO, several new infrared singularity operators appear; the one- loop soft gluon 
current S^{e,q 2 ) [12], 
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and the hard radiation functions H^e, g 2 ), which decompose into virtual contributions 



H,!-l(e, g 2 ) and real contributions H^ ; R (e, g 2 ). For the antenna functions presented here, 
the virtual contributions are given by, 
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while the real radiation contributions are, 
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Combining the real and virtual hard radiation functions yields the hard radiation terms 
from the two-loop virtual corrections [4,5,44,46,47], 



qq,N 2 



(2) 



1 



r(2) 



'/</: 



N 2 "gg.l/JV** 



+NN F H% NNF (e,f) + ^H^ F/JV ( e ,^) , 
H( 2 ) (e ,, 2 ) = ^h£> 2 ( e ,, 2 ) + JV^H£> M 2 ) + ^ H £> /Ar (e,, 2 



2A j_ ^>Tl( 2 ) 



19," F 

"■gg,N^ 
[(2) 



H^)( e ,g 2 ) = Ar 2 H^ 2 ( e , g 2 ) + AriV^^^^^^ 



+iV 2 Hi 2 ) ( £ ,, 2 ), 



(4.5) 



(4.6) 



(4.7) 
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HS^(^ 2 ) = ^K 2 l 4 ( e ,, 2 )+H^ (e , (? : 

- ^ /409 lbr 2 7 \ 2 2e 
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3Csl (9 2 )" 2e , (4.9) 



HSv„,v(c,« 2 ) = -HS> (e.g 2 ), (4.12) 
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The relation of these operators to the singularity structure of physical multi-parton two- 
loop matrix elements was discussed in [4,5,44,46,47] 

Finally, the NNLO singularity structures also contain the QCD /3-function (2.13) and 
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the collinear coefficient K. In a colour-ordered decomposition, these are 



and 



00 = b N + b ,FN F with 6 = ^ , &o,F = -7: 

6 3 



67 7T 2 

K = k N + k Q F N F with fc = — - — , fco,F 

18 6 



(4.21) 
(4.22) 



5. Quark-antiquark antennae 

The quark-antiquark antenna functions are derived by appropriately normalising the colour- 
ordered QCD real radiation corrections to 7* — ► qq, described to NNLO accuracy in [44]. 

The overall normalisation is given by defining the tree-level two-parton quark-antiquark 
antenna function 



1 



The one-loop two-parton quark-antiquark antenna is then: 
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with 



Voles (A\{s 12 )) = 2lg(e,s 12 ), 
^mte (^(si 2 )) = -4. 

5.1 Three-parton tree-level antenna functions 

The tree-level three-parton quark-antiquark antenna is: 



1 ^£13 £23 2 «12Sl23\ 



S123 V s 23 «13 S13S23 / 

yielding the integrated antenna function according to (2.11): 
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This antenna can be split symmetrically into two sub-antennae, which coincide with 
the q — > qg dipole functions: 

A§(1,3,2) =a°(l, 3,2)+ a°(2, 3,1) , (5.9) 

with 



4(1, 3, 2) = ^ + 2 g ^ + (9(e) . (5.10) 

«123 \Sl3 («13 + S 23 )Sl3/ 

5.2 Three-parton one-loop antenna functions 

At one loop, one finds three different three-particle antenna functions, corresponding to 
the leading and subleading colour structures A\{l q , 3 g , 2 q ), A\(l q , 3 g , 2 q ) and to the con- 
tribution from a closed quark loop A\(l q , 3 g , 2 q ). 
Introducing 

TT 2 

R{y,z) = log y log z - log y log (1 - y) - log 2; log (1 - z) + — - Li 2 (y) - Li 2 (z) (5.11) 
and yij = Sij/si2 3 , the one- loop antenna functions are given by: 
Voles (Al(l q ,3 g ,2 q )) = 2 (l«(e, si 3 ) + lg(e, s 23 ) - lg(e, a 123 )) A° 3 (l, 3, 2) , (5.12) 
finite (4(l„3 ff ,2 ? )) = - (fl(yi 3 , 2/23) + | log yi 3 + | logy 23 ) A§(1,3,2) 
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Voles (A&l q ,3 g ,2 9 j) = 2 (l« F (e, s 13 ) + I^ F (e, s 23 )) A°(l,3,2) , (5.16) 

•Fmi*e(A£(l„3 ff ,2 ? )) = 1 (logy 13 + logy 23 ) A§(1, 3, 2) . (5.17) 
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Note that application of the finite-operator in the above expression yields only the C(e )- 
terms of the antenna functions. These antenna functions contain higher powers in e as 
well, and these are relevant to the integrated antennae listed below. 
The integrated antennae are defined in (2.35). They read: 



Al(s 123 ) = (S123) 
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31 
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«4$(si23) = (-123)" 
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5 TT 2 
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"T + T + 7C3 



+ 



105 27^ 2 27, 7vr 4 ', 
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^3( s 123) = (-123)" 
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1 1 1 
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TV 
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with 



PoZe- (^3(5123)) 

Finite (^(s^)) 
Pok- (^(-m)) 
Finite ^3(^123)) 
PoZes (^(-123)) 
Finite (Al(e, S123)) 



-^( Sl23 ) (21^(6,^23) +^(S123)) + ( S 123)- e lS ) (e,Sl23) 

-Hg i ) A (e, Sl23 ) + S^ ) (e, Sl 23) , (5.21) 
6581 787vr 2 17tt 4 143 
48 + 96 + 360 + ~3~^ 3 ' 

-^(-123) (21^(6,^23) + A%(a 123 j) - H^(e, Sl23 ) , 

845 217vr 2 9tt 4 75 
~~8~ + 32 + lo~ + Y^ 3 ' 



2&o,. 



(S123) e liy(e,si23) 



e 

109 8 
~2T ~ 3^ 3 ' 



H« (-123) , 



(5.22) 
(5.23) 

(5.24) 
(5.25) 
(5.26) 



5.3 Four-parton tree-level antenna functions 

The tree-level four-parton quark-antiquark antenna contains three final states: quark- 
gluon-gluon-antiquark at leading and subleading colour, A® and A® and quark-antiquark- 
quark-antiquark for non-identical quark flavours B® as well as the identical-flavour-only 
contribution C4. The quark-antiquark-quark-antiquark final state with identical quark 
flavours is thus described by the sum of antennae for non-identical flavour and identical- 
flavour-only. The antennae for the qggq final state are: 



A° 4 (l q , 3 g , 4 g , 2g) = a°(l, 3, 4, 2) + a°(2, 4, 3, 1) 



(5.27) 



A° 4 (l q , 3 g , 4 g , 2g) = 52(1, 3, 4, 2) + a°(2, 4, 3, 1) + a° 4 (l, 4, 3, 2) + a°(2, 3, 4, 1) , (5.28) 
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where the sub-antennae are given by 
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[si2 - S34] + — + 0(e) 
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(5.30) 



In A® the gluonic emissions are colour-ordered, while in A® the gluons are photon-like, 
implying no ordering. Because of colour-ordering, A\ can be used with a single ordered 
phase space mapping. In contrast, A\ can not be used with a unique ordered phase space 
mapping. The above decomposition into a\ yields however ordered terms, since the com- 
bination a2(l, 3, 4, 2) + a2(2, 4, 3, 1) contains only single emission singularities in I/S13 and 
l/s24> corresponding to the ordered (1,3,4,2) phase space mapping. On the other hand 
52(1,4,3,2) +04(2,3,4,1) contains only single emission singularities in I/S14 and l/s23> 
corresponding to the ordered (1,4,3,2) phase space mapping. Since the decomposition of 
A® is symmetric, all four a\ yield identical integrals if integrated over the tripole phase 
space. It should be noted that it is not possible to analytically integrate an individual 
64 over the tripole phase space using the reduction and integration techniques described 
in [31], since the extra polynomial denominators present there enlarge the set of basis in- 
tegrals considerably. When the four a\ are added together these polynomial denominators 
cancel, and the tripole integrals can be carried out. 

The integrals of these antenna functions are according to (2.23): 
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with 



Voles (.42(51234)) = 2 1$ (e,Si234 



-2e" e7 



1 j + koj I^ ; (2e,si234j 
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— H^ A (e, S1234) — Sy'(e, S1234) , 
811037 23217T 2 13vr 4 4217 
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Finite (42(si234)) 
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(5.33) 
(5.34) 
(5.35) 
(5.36) 



32 16 10 
The non-identical quark antenna is: 

B° 4 (l q ,S q/ ,4g',2g) = 62(1, 3, 4, 2) +62(2, 3, 4,1) +62(1,4, 3, 2) +62(2, 4, 3,1) , (5.37) 

with a sub-antenna function given by 

1 f 1 
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In .B4, the secondary quark emission is ordered, such that a single ordered phase space 
mapping can be used. 

This subtraction term yields the integral 



#4(31234) = (S1234) 
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(5.39) 



with 



Voles (B2(si23 4 )) 
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r(l - e) 



- H R,BMl234) , 
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The identical-flavour-only quark-antiquark-quark-antiquark antenna is: 
Ct(l q , 3„ 4 g -, 2,-) = c2(l, 2, 3, 4) + cg(l, 4, 3, 2) , 



(5.40) 
(5.41) 

(5.42) 
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with 
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S13S| 4 S12S13 
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+ 7i " [-S12S14 - S12S34 - s\ 2 + S13S24I 

2S23S123«234 

1 r 2i S 13 

+- [S12SU + S12S34 + s l2 + S13S24 
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+ — 2 [Sl2«24 + S14S24] + 7^r—7— [~ s i2 + sw] + 0{e) } . (5.43) 



It integrates to 



C 4 ( S 1234) = j ('1234 ) 7 ri6 + Y"2 C3 J 



with 



, 339 17vr 2 21 2vr 4 \ 

+ (-^2- + ^r + T C3 -i5-' + 0(e) 



PoZes (C 4 °( Sl 234)) = -^Hg^e, Sl234 ) , (5.45) 

frM, 339 177r2 21 7T 4 , 

^mte (C 4 °( Sl2 34)) = + -^g- + yC3 - ^ • (5.46) 

All antenna functions listed in this section agree with the four-parton matrix elements 
in the Appendix of [35], taking account of the different normalisation used here. They 
also all agree with [Al]. 



(5.44) 



6. Quark-gluon antennae 

The quark-gluon antenna functions are obtained from the QCD real radiation corrections 
to the decay of a heavy neutralino into a massless gluino and a gluon, \ ~^ 99-, which is 
described in detail in [46]. 

The overall normalisation is given by defining the tree-level two-parton quark-gluon 
antenna function 

V° 2 (s 13 ) = 1 . (6.1) 

In this equation, and in all subsequent equations in the section, we label the primary quark 
momentum as (1) ? and the momenta of gluons or of a secondary quark-antiquark pair 
as (3)j, (4)j and (5)^. This non-consecutive labelling of momenta is introduced in view 
of applying the quark-gluon antenna functions in an actual calculation, where they will 
always appear in a pair: quark-gluon antenna and antiquark-gluon antenna, with (l) g and 
(2)q denoting the primary quark and antiquark momenta. 
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The one-loop two-parton quark-gluon antenna contains two contributions, correspond- 
ing to the different colour and flavour structures: 



Vl(s 13 ) = 2(s 13 ) 
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PoZes (^(si 3 )) 
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(6.2) 
(6.3) 

(6.4) 
(6.5) 
(6.6) 

(6.7) 



The pole terms in the above expression have to be compared to the pole terms of the 



(i) 

qq 



one-loop correction to the quark-antiquark antenna function A\ in (5.3), containing 21 
The factor 4 in (6. 4), (6. 6) appears since the tree level quark-gluon antenna function D^s^) 
contains two distinct quark-gluon antennae, in contrast to the single quark-antiquark an- 
tenna contained in A^is^), as will be seen below in constructing the three-parton tree- level 
antenna functions. 

6.1 Three-parton tree-level antenna functions 

The tree- level three-parton quark-gluon antenna contains two final states: quark- gluon- 
gluon and quark-quark-antiquark. The antenna corresponding to the first final state is: 
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Its integrated form is 

P^(S134) = 2 (S134)- 
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Voles {V° 3 (s lu )) = -4lW(e, Sl34 ) , 
Finite (T>%(s lu )) = ^ . 
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(6.10) 
(6.11) 
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This tree-level antenna function contains two antennae, corresponding to the configura- 
tions: (gluon (3 9 ) radiated between quark and gluon (4 g )) and (gluon (4 ff ) radiated be- 
tween quark and gluon (3 S )). The separation between these is not free from an ambiguity, 
since the collinear limit of the two gluons has be split between the two configurations. We 
decompose 

Dg(l, 3, 4) = 4(1, 3, 4) + 4(1, 4, 3) , (6.12) 
where the sub-antenna is given by 



dg(l, 3,4) = ^ + susu + sl + s_^ + 5 + . , + 

Si 3 4 \ S13S34 S13 «34 2 2 



The function could be further decomposed into two dipoles if needed for configurations 
discussed in Section 2.3.2. 

The tree-level three-parton quark-gluon antenna corresponding to the quark-quark- 
antiquark final state is: 

E$(l q ,3 q ,,4 q ,) = -L (il±A± + Sl3 + Su \ +0 (e). (6.14) 

s 134 V s 34 / 

Phase space integration yields: 
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with 



Voles (f§(si34)) = -4I^ F (e, s 134 ) , (6.16) 
finite («? 3 °(si34)) = -1 • (6.17) 

6.2 Three-parton one-loop antenna functions 

At one loop, the correction to the quark-gluon-gluon antenna contains a leading colour 
term D\(l q , 3 g , 4 g ) and a quark loop term D\(l q , 3 g , 4 ff ). These read: 
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+^1f(^ *34) - 21^ (e, ai34 ) D° 3 (l, 3, 4) , 



(6.20) 



Finite (b\{l q ,3 g A g )) = ^ (logyis + logj/u + 2 logy 3 4) 3, 4) - J- . (6.21) 
Integration of these antenna functions yields 
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Again, the factors of 2 in front of the infrared singularity operators and of the soft gluon 
current arise from the fact that the basic two-parton process contains two quark-gluon 
antennae. 

At one loop, the correction to the quark-quark-antiquark antenna contains a leading 
colour and a subleading colour term E\(l q , 3 q ', 4 q i), El(l q ,3 q ',4 q ') as well as a quark loop 
term E\(l q , 3 q >, 4 q i). These read: 
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Poles (^(l,^,^)) = 2 (lW(e, S3 4)) £§(1,3,4) , 
finite (^(l,^,^)) = -4E§(1,3,4) , 
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= -4I^ F (e, Sl3 4)^3 (l,3,4) , 
finite (El (l q , 3 q ,, 4,-0 ) = (-y + | log y 3 A £§(1, 3, 4) 

Integration of these antenna functions yields 
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The last contribution corresponds to the quark loop correction to the quark-quark-antiquark 
antenna, and enters into physical cross sections multiplied with a factor Np. In the case 
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of quark-gluon antenna functions, there is no contribution with this factor coming from 
four-parton tree-level antennae. Therefore, Voles{£\) contains the full infrared singularity 
structure of the two- loop two-parton quark-gluon antenna at Np. Note that the full hard 

(2) 

radiation factor ET^ n2 (and not only the virtual hard radiation) is present here. 
6.3 Four-parton tree-level antenna functions 

The tree-level four-parton quark-gluon antenna contains two final states: quark-gluon- 
gluon-gluon, and quark-quark-antiquark-gluon at leading and subleading colour, E® 
and E®. The antenna for the qggg final state is: 
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In Dl(l, 3, 4, 5), the gluonic emissions are colour-ordered. However, as explained in [46], 
both possible colour-orderings (normal and reverse) are contained in the same colour- 
ordered squared matrix element. Moreover, since this squared matrix element originates 
from an amplitude which is cyclic in the partonic colour indices, it also contains singular 
limits in the momentum configuration (5 g , l q , 3 g ), as outlined in detail in Section 8.3.2 be- 
low. Using this antenna function as NNLO subtraction term, this limit has to be accounted 
for properly. 

For the above reasons, it is not possible to use a single ordered phase space parametri- 
sation in the numerical implementation of D A (l, 3, 4, 5). Instead, D 4 (l, 3, 4, 5) has to be 
split into different ordered pieces by repeated partial fractioning on pairs of invariants. 
Since this decomposition is not symmetric, the individual pieces can no longer be inte- 
grated analytically. Therefore the full D 4 (l, 3, 4, 5) and not the sub-antennae must be used 
as a subtraction term. 

The integral of this antenna function is 
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The leading colour and subleading colour quark-quark-antiquark-gluon antennae are: 
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Concerning the phase space mapping, E® can be decomposed into ordered contributions, 
which can use an ordered parametrisation of the phase space, by repeated partial fractioning 
of the invariants (like in A®). Due to the lack of symmetry, the result is considerably larger 
than the above expression and not very instructive. It is therefore not quoted here. E® is 
already ordered (since quark (1) decouples from all singular limits). 

In the leading colour term, the gluonic emission is colour-ordered (either in the q, q'- 
antenna denoted here, or in the q, q'- antenna). The second colour-ordering is obtained by 
interchanging q' <-> q' . Note, in particular, that gluon emission inside the q', q'- antenna 
is only possible at subleading colour (and constitutes in fact the full subleading colour 
contribution) . 
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Integration yields 
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7. Gluon-gluon antennae 

The gluon-gluon antenna functions are obtained from the QCD real radiation corrections 
to the decay of a massive Higgs boson into two gluons, H — > gg, which is described in 
detail in [47]. 

The overall normalisation is given by defining the tree-level two-parton gluon-gluon 
antenna function 

^(si 2 ) = 1 • (7.1) 

The one- loop two-parton gluon-gluon antenna contains two contributions, correspond- 
ing to the different colour and flavour structures: 
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The pole terms in the above expression have to be compared to the pole terms of the 
one-loop correction to the quark-antiquark antenna function A\ in (5.3), and the quark- 
gluon antenna functions T>\ in (6.4) and T)\ in (6.6). The factor 4 in (7. 4), (7. 6) appears 
since the tree-level gluon-gluon antenna function contains two distinct gluon-gluon 
antennae. This situation is like for V>\, which contains two quark-gluon antennae, but in 
contrast to the single quark-antiquark antenna contained in A®- 



7.1 Three-parton tree-level antenna functions 

The tree- level three-parton gluon-gluon antenna contains two final states: gluon-gluon- 
gluon and gluon-quark-antiquark. The antenna corresponding to the first final state is: 
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It yields the integral: 
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As can be seen from the pole structure, this tree-level antenna function contains three 
antenna configurations, corresponding to the three possible configurations of emitting a 
gluon in between a gluon pair. The separation between these is not free from an ambiguity, 
but is in fact fixed by the decomposition used in the case of the quark-gluon-gluon antenna 
in (6.13) above. We decompose 
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The tree-level three-parton gluon-gluon antenna corresponding to the gluon-quark- 
antiquark final state is: 
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Its integrated form reads: 
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7.2 Three-parton one-loop antenna functions 
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At one loop, the correction to the gluon-gluon- gluon antenna contains a leading colour 
term ^3(51, 52, 53) and a quark loop term ^3(51, 52, 53)- These read: 
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Integration of these antenna functions yields 
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with 
Voles (^(si2s)) 
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The factors of 3 in front of the infrared singularity operators and of the soft gluon current 
arise from the fact that the tree-level three-parton antenna function contains three antenna 
configurations. 

At one loop, the correction to the quark-quark-antiquark antenna contains a leading 
colour and a subleading colour term G\(l g , 3 q , A q ), G\(l g , 3 q , A q ) as well as a quark loop 
term Gl(l g ,3 q ,A q ). These read: 



Voles {Gl(l g , 3„ 4,-)) = 2 (l«(e, s 13 ) + I«(e, s 14 ) - 2 I«(e, ai34 )) G§(1, 3, 4),(7.28) 

/ 5 5 

Finite (Gl(l g ,3 q ,4q)) = -( R(y 13 , y u ) + R(yu, 2/34) + ^ log 2/13 + ^ log yu 

+ ^^2/34-^^(1,3,4)-^— 

Voles (Gi(l 9 ,3 9 ,4 9 -)) = 2 (lg(e,s 3 4)) G°(l,3,4) , 
Finite (Gl{l g ,3 q ,4g)) 
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9 +|logj/34 + ^logyi 3 + ^logyi4)G3(l,3. I.). (7.;!:!) 
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Integration of these antenna functions yields 
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In contrast to the Np one-loop three-parton quark-gluon antenna function £3, which had 
no corresponding tree-level four-parton antenna function with the same colour and flavour 
structure, Q\ is complemented by an antenna function containing two quark-antiquark 
pairs, H®, which is derived below. Correspondingly, Q\ does not contain the full iVj,-terms 
for gluon-gluon final states. 

7.3 Four-parton tree-level antenna functions 

The tree- level four-parton quark-gluon antenna contains three final states: gluon-gluon- 
gluon-gluon, F®, and gluon-quark-antiquark-gluon at leading and sub leading colour, G 4 
and (5° and quark-antiquark-quark-antiquark, H®. The antenna for the gggg final state is: 



^4(91,92,93,94) 



/ 4 °(1, 2, 3, 4) + / 4 °(4, 3, 2, 1) + / 4 °(2, 3, 4, 1) + / 4 °(1, 4, 3, 2) 
+/ 4 °(3, 4, 1, 2) + / 4 °(2, 1, 4, 3) + / 4 °(4, 1, 2, 3) + / 4 °(3, 2, 1,4) , (7.43) 
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with 
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In F®, the gluonic emissions are colour-ordered. Since the original colour structure is a 
trace over the gluon colour indices [47], is symmetric under cyclic interchanges of the 
momenta. Therefore, each pair of adjacent momenta can act as hard emitter pair for the 
antenna function. F® thus contains four different colour-ordered antennae. For numerical 
implementation, these have to be separated from each other by repeated partial fractioning 
of the associated invariants. This fractioning is not made explicit in the above expression 
fl 

The integral of this antenna function is 
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with 
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162 27 360 6 

The leading colour and sublcading colour gluon-quark-antiquark-gluon antennae are: 



G° 4 (l g ,3 q ,Ag,2 g ) = ^(1,3,4,2) + 5 °(2,4,3,1) 



(7.48) 



G° A (l g , 3„ 4,-, 2 fl ) = 5°(1, 3, 4, 2) + ~g° 4 (l, 4, 3, 2) + g° 4 (2, 3, 4, 1) + g° 4 (2, 4, 3, 1) , (7.49) 
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In the leading colour contribution G4, the gluonic emissions are colour-ordered in between 
the quark-antiquark pair. However, like all gluon-gluon antenna functions at leading colour, 
G4 is cyclic in the colour indices, such that it also contains a configuration where the two 
gluons form the hard emitter pair, emitting the quark-antiquark pair inside the antenna. 
As before, the different antenna configurations can be separated from each other for the 
numerical integration by repeated partial fractioning of the invariants. 
Integration yields 



<? 4 (si234) = 77 (S1234) 



-2c 



155 1 / 523 79vr' 



+ - 



12 + 36 



2e 3 18e 2 e 

16579 13857T 2 



+ - 



81 



+ 



108 



+ 37C 3 + 0(e) 



, (7.52) 



<?4(' s 1234) = (S1234) 



-2c 



1 



41 



+ - 



1 / 1327 vr 2 



3e 3 18e 2 e 



108 + 2 



- 56 - 



with 



+ 



4864 41vr 2 86 , . 



Voles (g° 4 (s 12U )) = 2 [21 W (e, 



gg v-> s 1234, 
-eyr(l-2e) 



r(i 



2lW F (e, S i234) 



2i« (2 e , si234 ; 



finite (£4(51234)) 

PoZeS (^(51234)) 

Finite ^°(si234 



16579 155tt 2 

+ 



162 



48 



649 

"36~ C3 ' 



-H^(e, ai2 34), 
4864 697tt 2 



81 



+ 



85 

1 C3 • 

216 9 s 



The quark-antiquark-quark-antiquark antenna is 
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where only different quark flavours need to be considered, since the identical flavour con- 
tribution to this final state is finite. 

The integral of this antenna term is 
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The contribution from identical quarks only to this final state is finite, and thus no gluon- 
gluon antenna function is defined for two identical quark-antiquark 



pairs. 
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8. Infrared limits of the antenna subtraction terms 



The antenna subtraction terms denned in the three previous sections encapsulate all single 
and double unresolved limits of tree-level QCD matrix elements and all single unresolved 
limits of one-loop QCD matrix elements. In this section, we list the behaviour of all 
tree-level three-parton and four-parton antenna functions. Using this information, it is 
then possible to employ these antenna functions in the construction of complete infrared 
subtraction functions for QCD matrix elements at NNLO. 

8.1 Generalised collinear and soft factors 

The factorisation properties of tree-level QCD squared matrix elements at NLO and NNLO 
have been investigated in detail in [17-20,38,40]. At NLO, only a single particle can 
become unresolved, either soft or collinear. In these limits, the (m + l)-parton matrix 
element factorises into a reduced m-parton matrix element times a soft eikonal factor or 
a collinear splitting function. At NNLO, two particles can become unresolved in several 
possible configurations: double soft, soft/collinear, double single collinear, triple collinear. 
In each of these limits, the (m + 2)-parton matrix element factorises into a reduced m- 
parton matrix element times a generalised double unresolved factor (double soft factor, 
soft/collinear splitting function, double single collinear splitting function, triple collinear 
splitting function). In the following, we list all generalised single and double unresolved 
factors. 

8.1.1 Single unresolved factors 

When deriving limits of the three-parton antenna functions in single unresolved configura- 
tions, we encounter well-known soft eikonal factors when a gluon is soft and three different 
Altarelli-Parisi splitting functions when two partons are collinear. These are listed below. 

When a soft gluon (b) is emitted between two hard partons (a and c), the eikonal factor 
S a b c factorises off the squared matrix element: 

2s ac 

Sabc = • (8.1) 

SabSbc 

When two partons become collinear, we have different splitting functions corresponding 
to various final state configurations: a quark splits into a quark and a gluon (P qg ^Q), a 
gluon splits into a quark-antiquark pair (P qq _,Q) or a gluon splits into two gluons (Pgg^c)- 
These are given by: 

'1 + (1- zf -ez 2 " 



Pqg^QW 



Z 

z 2 + {l-zf 
1 -e 



P 99 ^G{z) = 2 {^- z + i-I + *(1 - Z )j . (8.2) 

In these equations, z is the momentum fraction of one of the collinear partons and the label 
q appearing in these splitting function can stand for a quark or antiquark: P qg -,Q = P qg ^Q 
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by charge conjugation. Q or G appearing in the collinear splitting functions denotes the 
parent particle of the two collinear partons i and j. In the discussion below, the momentum 
associated with this parent particle will be denoted by (ij). 

8.1.2 Double unresolved factors 

To describe the limiting behaviours of the four-parton antenna functions in double unre- 
solved configurations, we require generalised soft and collinear factors. These were first 
derived in [17-20], and are listed below. We follow largely the notation of [18]. 

1. Double soft factors 

When two colour connected gluons (b) and (c) are simultaneously soft, the double 
soft gluon function for four partons in the final state is given by 

2s ld 

Sabcd(.Sad, Sab, Scd, Sbc, Saba Sbcd) — 



S abSbcdS abcS cd 

| 2s ad ( 1 + 1 + 1 _ 4 \ + 2(1 -e) / + £cd__{ 

Sbc \ s ab^cd S ab S bcd S cd S abc S abc S bcd J S^ c \Sabc S bcd 

(8.3) 

Here a and d are the hard partons surrounding the soft pair. 

In the case when two unconnected gluons become simultaneously soft the correspond- 
ing double soft gluon factor is the product of two eikonal factors S abc , given in (8.1). 
For the emission of a soft quark-antiquark pair resulting from the splitting of an 
intermediate gluon emitted itself from a primary pair of hard partons, the soft factor 
takes the form, 

2 

S<xb{Cqi dq) = - s -. ; r~. ; r (s ab S cd S ac S bd S bc S ad ) 
Si d {s ac + S ad ){Sbc + S bd ) 

+ 4- ( -^±- 2 + - f SbcSb \ 2 ) , (8.4) 

s z cd V («ac + s ad y {s bc + s bd y J 

with a and b the adjacent hard partons. Note that this formula is different from the 
one that can be found in [19] for this limit, where the second line of the above formula 
is absent. 

2. Soft-collinear factors 

When a gluon (a) is soft, partons (d,a,b,c) are colour connected and partons b and c 
are collinear then the limiting behaviour of the squared matrix element is described 
by the soft-collinear factor, which is subsequently multiplied with the appropriate 
simple collinear splitting function 

c ( \ ( s bd + s cd ) ( s ab + zs bc \ 

bd;abc{Z, S a b, Sbc, S abc , S ad , S bd , S cd ) — I Z H I . (S.t>J 

SabSad \ S abc J 

In here, z is always the collinear momentum fraction of parton b in the collinear pair 
(6c), with b being colour-connected to the soft parton a. 
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3. Triple collinear splitting functions 

There are five different triple collinear splitting functions depending on the nature of 
the partons which become collinear. 

The colour-ordered triple collinear splitting function P ggg ^c is given by 

Pabc^G(w,x,y,s ab ,s bc ,s abc ) = 2 x I 

(1 - e) {xs abc - (1 - y)s bc ) 2 , 2(l-e)s 6c , 3(1 - e) 



+ s lb s lbc C 1 " V) 2 + s abS 2 abc + 2s 2 



abc 



, ' (1 ~ yO- - y)) 2 2 x 2 + xy + y 2 | xw-x 2 y-2 | 



sabSabc V y w ( l ~ w ) ~ 1 ~y y(i - y) (i - y) 

+ 1 / 3x 2 _ 2(2 - w + w 2 )(x 2 + w(l - w)) 1 



2s afe s fec V y(i-y) (i-yX 1 -^), 

+ (Safe <-> Sbc,^ <-> y) , (8.6) 

with iu, x and y being the momentum fractions carried by the collinear gluons a, b, c, 
and with w = (1 — x — y). This splitting function is symmetric under the exchange 
of the outer gluons (a and c), and contains poles only in s ab and s bc . 
Pqg ig2 -^Q is the triple collinear splitting function obtained when quark q is collinear 
to two gluons, g\ and gi , which are not colour connected, i.e. behave like photons [17]. 
It reads, 

Pqgig2^>Q{ w i x i y-i s qgi > s qg2 j s qgig2) = 

+ — (l + w 2 - e(x 2 + xy + y 2 ) - e 2 xy) 

qgi s qg2 x y 

H — (w(l - x + e 2 xy) + (1 - y) 3 - e(l - y)(x 2 + xy + y 2 ) + e 2 xy) 

s qgi s qgig2 x y 

- %^ f (i - ^ - e ) + ~ s ^ x ~ y)> ( 8 - 7 ) 

with x,y and w the momentum fractions carried by the collinear particles and with 
w = (1 — x — y). For the case when the two gluons are instead colour connected, the 
colour-ordered splitting function reads, 

Pqgig2^>Q{ w i x i y-i s qgi > s qg2 > s 9i92 j s qgig2) = 

, 1 + w 2 l + (l-y) 2 \ /to 1-y 



+ — 



y (1 — w) J \y 1 — w 

[1 - y) 3 + u>(l — x) — 2y 



+ 



y(l - it>) 
' 2(1 -y)(y-w) \ 2 \ 

1 - y) 2 (2 - y) + x 3 + 2xw -2-y 
y(i - w) 
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(zw - 2/ - 2yw) \ 
y(l -w) J 

:(-,_,( 2(xs qgig2 - {l-w)s qgi ) 2 
'[ s 2 s 2 (1 - w) 2 

+J— Limn. + (1 _ e) w + (3 _ e) U (88) 

Similarly the clustering of a gluon with a quark-antiquark pair into a parent gluon 
again has two distinct functions depending on whether the final state is colour-ordered 
or not, In the colour connected case, the gluon is emitted outside the quark-antiquark 
pair, and one obtains the colour-ordered splitting function, 



^gqq— >g(^> Vi Sgq-i ^qqi Sgqq) 



S f y _L (1 c \ Sqc l _l_ d A\ ^ {xSggg (1 W)s gq ) 

s 2 V~ ( 'V- 1 ' s 2 s 2 (I - w) 2 

b gqq \ b qq b gq / b qq b gqq\ L w ) 

i ( (i-y) 9 2x(i-j/)(j/- w ) 

+ — n V -y-2w-e - — — 

SgqSgqq \W(1 - W) (1 - €)w{l - W) 

' x ((1 - w f - w 3 ) 2x 2 (1 - yw - (1 - y)(l - w)) 
w(l — w) (1 — e)w(l — w) 



2 



+ 



+ 1 / (1 + + 4:xw) + 2x (w(a- ;/) - y(l + w)) \ 

SqqSgqq\ w(l - w) (1 - e)w(l - w) /' 

while in the case where the gluon is emitted between the quark-antiquark pair, which 
is subleading in colour, one obtains a QED-like splitting function, 

Pqgq^ciw i x iVi s qgi s gqi s qqi s qgq) = 

-^(<l- e) £* + l) + ^- ((! + *>) 

b qgq \ b qg J b gq b qg \ L fc 

— ( 1 + 2x + e - y J + (s 99 <-> s M , w <-> y). 

s qg s qgq V v 1 e J / 



(8.10) 



Lastly, the clustering of a quark-antiquark pair (q'q r ) and a quark (q) to form a parent 
quark Q with the same flavour as q. The splitting function depends upon whether or 
not the quarks are of identical flavour, 

D _ pnon— ident. "<?<?' pident. /o i i \ 

Uq'q'^Q ~ ^qq'q'^Q ~ JV^Q' ^-^J 

where = 1 for identical flavour quarks. If quarks q±, q' 3 and q' A are clustered to 
form Q, 

Pqq'qiL^Q ( w i x i U i s qq' i s q'q' i s qq'q') 

2s q g> \ 2 (xSq q ,q f - (1 - w)s q -q l ) 2 



S 2 ,_, V ^ ^ ^ s a'a' J S 2 ,-,S 2 ,_,(1 — W) 2 

qq q x 9? ' g'q qq 'q \ ' 

+ A + *» + (* + u,)° - e(l - • (8.12) 

Sq'q'Sqq'q' V I 1 ~ / 
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When the flavours of the clustering quarks are the same, there is an additional contri- 
bution coming from the interference terms of the four-quark matrix elements, which 
reads 

P qfqL Q K X, y, S qljl , S qq > , S qq , q , ) = 

_ (l-€) /2^ |2| \ 1 / x(l + x 2 ) , / 2(l-y) |1|c 

S lq'q' \ S l'f / 2s qq' S q'q' \( l ~ V){ 1 ~ W ) ^ V i 1 ~ W ) 6 

1 fl + x 2 2x f(l~w) 2 M 2x ^ 
+ 7: r + 71 T " 6 Ti V + (1 + x) + r + e l - w) 

Sq>g>S gq <q' \(l-y) (1 - W) \ (1 - y) (1 - W) 

+ (Sqq> ^ Sq'q',y w). (8.13) 

This identical flavour splitting function contains poles when q' clusters with both q' 
and q. It is symmetric under q\ <-> q' 3 . 

4. Double collinear splitting functions 

If two distinct pairs of partons become simultaneously collinear, the (m + 2)-parton 
squared matrix element factorises into the product of two simple collinear splitting 
functions (8.2) with the m-parton squared matrix element. 

8.2 Quark- ant iquark antennae 

The three-parton and four-parton quark-antiquark antenna functions were derived in Sec- 
tion 5 from the real radiation corrections to 7* — ► qq. Their behaviour in all limits where 
one or two partons become unresolved are described below. 

8.2.1 Three-parton antenna functions 

We have only one three-parton quark-antiquark antenna function: A®(l q ,3 g , 2 q ), which has 
the following limits: 

1. Soft limit: 

Ajj(l,3,2) 3 ^°Si 32 . (8.14) 

2. Collinear limit: 

4(1,3,2) ^ ^-P qg ^ Q (z). (8.15) 

■513 

where the momenta of the quark (l q ) and the antiquark (2 q ) can be interchanged. 

8.2.2 Four-parton antenna functions 

The NNLO real radiation corrections to 7* — ► qq yield four different four-parton quark- 
antiquark antenna functions: the leading and subleading colour quark-antiquark-gluon- 
gluon antennae A®(l q , 3 g , 4 S , 2 q ), A®(l q , 3 g , 4 S , 2 q ), as well as the antennae with two quark- 
antiquark pairs of different and identical flavour B®(l q ,3 q > ,4 q /,2 q ) and C4 (l q , 3 q , 4 q , 2 q ). 
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Their single and double unresolved limits are described in the following. We shall always 
restrict ourselves to the non-vanishing limits only. 

In ^2(1,3, 4, 2), the colour-ordering of the gluonic emissions ensures that only neigh- 
bouring partons can form singular configurations. While singularities are present in I/S13 
and I/S24, no singularities in 1/su and I/S23 appear. 

For this colour-ordered antenna function, we find the following non-vanishing double 
unresolved limits: 

1. Double soft and soft-collinear limits: 

4(M,4,2) 3fl ^T° S 1U2 , 

4(M,4,2) 39 ^ 112 *- S 1;3 42(*) — P qg -> Q (z) , 

■524 

4(M,4,2) ^-^l" 1 ' S 2 .M*) — Vo(^) • (8-16) 

«13 



2. Triple collinear limits 



^(1,3,4,2) lql H A9 P 13 ^ Q (w,z,y) , (8.17) 



with quark (l q ) and antiquark (2q) being interchangeable. 
3. Double collinear limits: 

4(M,4,2) l4 ^ hi J- P qg ^ Q {z) J- P^ Q (y) . (8.18) 

We find the following non-vanishing single unresolved limits: 

1. Soft limits: 

A° 4 (l,i,j,2) ^ S Uj A°(l,j,2), 

A° 4 (l,i,j,2) j ^° S ij2 A° 3 (l,i,2). (8.19) 

2. Collinear limits (quark-gluon) : 

4(1,3,4,2) * — P qg ^ Q (z) A°((13),4,2) , 

^(1,3,4,2) * — P qg ^ Q (z) ^(1,3,(24)) , (8.20) 

where the parent parton of the collinear partons (1) and (3) (or (2) and (4)) is denoted 
(13) (or (24)) in the three-parton antenna function A®. 

3. Collinear limit (gluon-gluon) : 

4(1,3,4,2) ^ — P gg ^G(z) 4(1, (34), 2) + ang. , (8.21) 

«34 

where (ang.) means that angular terms are also present here. These do however 
cancel after integration over the unresolved phase space. A detailed discussion of 
these angular terms can be found in Section 8.5 below. 
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The limits of the other colour-ordered subtraction term, ^(l, 4, 3, 2), can be inferred 
from the limits of ^(l, 3, 4, 2) listed above by interchanging the momenta of the two final 
state gluons (3 S ) and (4 S ). 

The subleading colour contribution ^2(1,3,4,2) is not colour-ordered, and therefore 
symmetric under the interchange of the two gluon momenta 3 and 4. We find the following 
non- vanishing double unresolved limits: 

1. Double soft and soft-collinear limits: 

I°(l,3,4,2) 3 ^lT° S 132 S U2 , 

j,2) 9i ^ h2 S 1;ij2 (z) — P qg ^ Q (z) , (8.22) 

Sj 2 

with i,j standing each for one of the two intermediate gluons. Quark (l g ) and anti- 
quark (2q) are interchangeable. 

2. Triple collinear limits: 

10(1, 3 , 4, 2) lql H 49 P lu ^ Q (w, x, y) , (8.23) 
with quark (l q ) and antiquark (2^) being interchangeable. 

3. Double collinear limits: 

Al(l,i,j,2) J- P qg ^ Q (z) J- P qg ^ Q (y) , (8.24) 

Sli S 2 j 

with i and j each standing for one of the two gluons (3 9 ) or (4 S ); quark (l g ) and 
antiquark (2^) can also be exchanged. 

The non-vanishing single unresolved limits are: 

1. Soft limits: 

i°(l,i,j,2) ^ S li2 A° 3 (l,j,2) , (8.25) 
where % and j represent the two gluons (3 g ) and (4 g ). 

2. Collinear limits: 

A° 4 (l,i,j,2) ^ — P qg -+Q(z) A° 3 ((U),j,2) . (8.26) 

Again, (l q ) and (2 q ) can be interchanged and % can stand for one of the two gluons (3 S ) or 
(4 S ). 

For the antenna function containing two quark-antiquark pairs of non-identical flavours, 
-B4(lg>3g',4q/,2q), we find the following double unresolved limits: 

1. Double soft limit: 

n 3„; -^0,4^,^0 

B 4 °(l,3,4,2) q -A S 12 (3,4). (8.27) 
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2. Triple collinear limit: 



$2(1,3,4,2) M M V P~ ^( X ,y) 
The quark (l g ) and antiquark (2^) momenta can be interchanged. 
54(1,3,4,2) contains only one single unresolved limit (3 ? || A q ): 

52(1,3,4,2) ^ — P qq ^ G {z) A§(1,(34),2) + ang. , (8.29) 

where (ang.) means that angular terms are also obtained here. 

Finally we consider the identical-flavour-only antenna function denoted by C®(l q , 3 q , A q , 2 q ). 
This function is proportional has only one non-vanishing double unresolved limit: triple 
collinear (2 q \\ 4 q || 3 ? ), 

Cj (1, 3, 4, 2) 2? ^t 4? \ Pgfyiw, x, y). (8.30) 

This antenna function has no further non- vanishing (double or single unresolved) limits, 
although terms in I/S24 and I/S34 are present in the function. In the respective limits, the 
coefficients of these terms vanish. 

8.3 Quark-gluon antennae 

The quark-gluon antenna functions were derived in Section 6 above. We describe their 
behaviour in all single and double unresolved limits in the following. Again, we will only 
list the non-vanishing contributions. In all these functions, the quark (l q ) can also represent 
an antiquark. 

8.3.1 Three-parton antenna functions 

The tree-level three-parton quark-gluon antenna functions represent two final states: quark- 
gluon-gluon (D®(l q , 3 ff , 4 ff )) and quark-quark-antiquark (E® (l q , 3 q >, 4^/)). 

The tree- level antenna function D®(l q , 3 9 , 4 ff ) contains two colour orderings correspond- 
ing to the following configurations: gluon (3 S ) radiated between quark (l g ) and gluon (4 ff ) 
denoted by dg(l,3,4) and gluon (3 S ) radiated between quark (l q ) and gluon (4 ff ) denoted 
by d§(l,4, 3), such that 

£>§(1,3,4)=^(1,3,4)+^(1,4,3) . 

The separation between these is not free from an ambiguity. It is constructed in such a way 
that the collinear limit of the two gluons has to be split between the two configurations. 
The decomposition used here is stated in (6.13). 

The simple unresolved limits of D®(1,3, 4) and ^(1,3,4) are stated below. For the 
sake of clarity, we also list some non-trivial vanishing limits: 



(8.28) 
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1. Soft limits: 

4(l,i,j) 3 ^0. (8.31) 



^ -*(!-,)), 



(8.32) 

The antenna function E®(l q , 3 9 ',4g/) has only one singular behaviour, when the quark- 
antiquark pair becomes collinear: 

E$(l, 3, 4) — P qq -+ G (z) . (8.33) 

■S34 

8.3.2 Four-parton antenna functions 

The NNLO real radiation corrections to x ~^ 99 yield three different four-parton quark- 
gluon antenna functions: the quark-gluon-gluon-gluon antenna contains only a leading 
colour (colour-ordered) term D^(l q , 3 g , 4 g , 5 g ), while the quark-quark-antiquark-gluon an- 
tennae have leading colour and subleading colour contributions E®{l q , 3 q /, 4 q i, 5 g ) and 
^4(lg)3 g ',4g/,5 g ). By construction, no identical-flavour contribution appears. The non- 
vanishing single and double unresolved limits of these antenna functions are described in 
the following. As before, some of the single unresolved limits also contain angular terms, 
which will be discussed in more detail in Section 8.5 below. 

The colour-ordered quark-gluon-gluon-gluon antenna function D®(l q , 3 g , 4 g , 5 g ) con- 
tains several different antenna configurations: due to the cyclic nature of the colour indices 
of the x ~^ 9999 matrix element [46], each pair of two neighbouring partons can represent 
the two hard partons forming the antenna. This behaviour is in contrast to the quark- 
antiquark antennae described in the previous subsection, where the colour-ordering en- 
sures that the primary quark-antiquark pair always forms the hard partons. Consequently, 



2. Collinear limits: 

D 3 °(l,3,4) 3 ^ —P g9 ^ G (z), 

d§(1.3,4) ^ J- (p^CO - 
d§(M,3)^ J- (p^C*) - 

S34 v 

D° 3 (l,i,j) ^ ^-P qg ^ Q (z), 

Slj 

djj(M,j) ^0. 
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Z?4 contains considerably more double and single unresolved limits than the four-parton 
quark-antiquark antenna functions studied above. 
For the double unresolved limits we find: 



1. Double soft and soft-collinear limits: 



D 4 °(l,3,4,5) 3 *^lr° <? 1345 , 
£ 4 (1,3,4,5) 4 *^T° S 1543 , 



£0(1,3,4,5) 



3 9 -+0,5 g ->0 



<Sl34 5*154 , 



(8.34) 



£0(1,3,4,5) ^ •" S^(z 

4 o ||5 o ,3 o ^0 



£0(1,3,4,5) 
£0(1,3,4,5) ^is-lT 



,0m o a e\ ^PgAg^O 



£^(1,3,4,5) 
£^(1,3,4,5) 



la 3o i So ^0 



£0(1,3,4,5) 39 " 4 ^° 



5l;345(^ 
S l 5;43l(2 
S , 3 ;45l(2 ; 
5'4;513(^ 
S'l;543(^ 



— Pqg^Q^-z) , 

815 
S45 

— P Q9^Q( Z ) 
813 

1 



qg- 



)(*), 



Sl5 

Sl3 
1 

S34 



(8.35) 



2. Triple collinear limits: 



■)0/l o A r\ 1 9ll 3 9ll 4 s 



£^(1,3,4,5) ^ 9 Pi34^oKx,y) 
£0(1,3,4,5) 19 P 154 ^ Q (u;,x,y) 



£^(1,3,4,5) 



135- 



£0(1,3,4,5) 39 ^4 5h P 



■ M :^ G {w,x,y) . (8.36) 

In these, the presence of the triple collinear limit (l q \\ 3 g \\ 5 g ) is particularly 
noteworthy. This limit appears due to the cyclic nature of the £4 antenna function. 
In contrast to the other two triple collinear limits involving (1 ? ), where the leading 
colour triple collinear splitting function Puj^q appears, this limit is controlled by the 
subleading colour splitting function -Pi35^Q, since the gluons (3 9 ) and (5 9 ) are not 
directly colour-connected. Applying £4 as antenna subtraction term to a physical 
multi-parton matrix element, special care has to be taken about this particular limit, 
which is a priori oversubtracted. 



3. Double collinear limits: 



£0(1,3,4,5) 



l 9 l|3 s ,4 9 ||5 fl 



£^(1,3,4,5) l9l|5 £l^ l|4s 



S13S45 
1 

S15S34 



Pqg-^Q(Z) Pgg^ G (y) 

p qg^Q.( z ) p gg^c(y) 



(8.37) 
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For the single unresolved limits we obtain: 
1. Soft limits: 



D 4 °(l,3,4,5) 3 -^°5i34 £>§(1,4,5), 
£>°(1,3,4,5) 4 -^ 5 3 45 ^(1,3,5), 

D° 4 (l, 3, 4, 5) 5 -^° S 154 D° 3 (l, 3, 4) . (8.38) 



2. Collinear limits 



£.0(1,3,4,5) * —P qg -> Q (z) £>§((13),4,5) , 

•513 

D2(l,3,4,5) * — Vq(2) £> 3 °((15),3,4) , 

^2(1,3,4,5) ^ — P^ G (z) D 3 °((l,(34),5) + ang. , 

£> 4 °(1,3,4,5) * — P^ G (z) J D 3 °((l,3,(45))+ang. . (8.39) 

The leading colour quark-quark-antiquark-gluon antenna function E®(l q , 3 q /, 4q>, 5 g ) 
is colour-ordered: the gluon 5 g is emitted only between the primary quark (l q ) and the 
secondary antiquark (4^/). The second colour-ordering E® (l q , 4^/, 3 q /, 5 ff ) is obtained by 
exchanging the secondary quark and antiquark momenta. 

For the double unresolved limits of this antenna function we find: 

1. Double soft and soft-collinear limits: 

£2(1,3,4,5) V ~^r"°S 15 (3,4) , (8.40) 

£0(1,3,4,5) 3 *' I|4 I^° 5 (3) J_ P 9 ^ G (z) . (8.41) 

S34 



2. Triple collinear limits: 



££(1,3,4,5) 1?II M V p^ dcnt -(^^,y) , 

£ 4 °(1, 3, 4, 5) 3 *'^J 59 P 543 _ G (™, x, y) . (8.42) 



3. Double collinear limit: 



E° 4 (l, 3, 4, 5) M5 ^' IIV _L_ P w -q(z) P q ^ G (y) ■ (8-43) 

S34S15 



For the single unresolved limits we obtain: 
1. Soft limit: 



P°(l, 3, 4, 5) ^ 5 154 E° 3 (l, 3, 4) , (8.44) 
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2. Collinear limits: 

£#(1,3,4,5) ^ — P q9 ^a(z) £>g(l,(34),5) + ang. , 

£2(1,3,4,5) ^ — P qg ^ Q (z) E$(l, 3,(45)) , 

£2(1,3,4,5) * — P^ Q (z)£ 3 °((15),3,4) . (8.45) 

In the subleading colour antenna E®(l q , 3 q i,A q i, 5 ff ) the gluon (5 ff ) is radiated between 
the secondary quark-antiquark pair(3 9 ', 4g/). As a consequence, no singular structure in- 
volves the primary quark (l q ), thus limiting the number of singular configurations con- 
tained in this antenna. It contains only one double unresolved configuration, triple collinear 

II 4q, || 5 g ): 

E° 4 (l, 3, 4, 5) 3 ^ |5fl P 5U _ G (w, x, y). (8.46) 

The single unresolved limits always involve the gluon (5 fl ). They read: 

1. Soft limit: 

E° 4 (l, 3, 4, 5) S 354 E° 3 (l, 3, 4) . (8.47) 

2. Collinear limits: 

£ 4 °(M,j,5) ^ X P qg ^ Q ( z ) £$(l,(i5), j) , (8.48) 
where i and j can both play the role of the quark (3 q /) or the antiquark (4^/). 
8.4 Gluon-gluon antennae 

The gluon-gluon antenna functions were derived in Section 7 above. Their behaviour in 
all single and double unresolved limits is summarised in the following, where only the 
non-vanishing contributions are given. 

8.4.1 Three-parton antenna functions 

There are two tree-level three-parton gluon-gluon antenna functions: gluon-gluon-gluon 
F®(l g , 2 g , 3 g ) and gluon-quark-antiquark G®(l g , 3 q , 4 q ). 

The tree-level antenna function F 3 (1 9 , 2 9 , 3 9 ) contains three antenna configurations 
since each pair of gluons can represent the hard partons, emitting the remaining third 
gluon. Each of these configurations is denoted by f®{i,j, k): 

F 3 °(l, 2, 3) = / 3 °(1, 2, 3) + / 3 °(1, 3, 2) + / 3 °(2, 1, 3) . 

As with the quark-gluon-gluon antenna function, this decomposition into antenna config- 
urations is not unambiguous, the decomposition used here is stated in (7.13). 
The simple unresolved limits of (1, 2, 3) and / 3 (1, 2, 3) are: 
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1. Soft limits: 



F 3 °(l,2,3) ^°5 213 , 
F 3 (l,2,3) 2 -^> 5 12 3, 
F 3 °(1,2,3) 3 -^°S 132 , 
/ 3 °(1,2,3) ^°S 123l 



/ 3 °(1,2,3) ^O, 

/ 3 °(1,2,3) ^°0, (8.49) 
where we also include the non-trivial vanishing limits for clarity. 
2. Collinear limits: 

F 3 (l,2,3)*— P^Z), 



Sl2 

^(1,2,3)^1^) 

«13 



-z(l-z) 



1 - z 



F 3 °(l,2,3)^ -=-J^ G (z), 
^(1,2,3) ^ J- (p^ g (z) 
/ 3 (1,2,3) — > — [P gg ^ G (z) z(l -z)\ , 

•523 V Z J 

/ 3 °(1,2,3) ^O. (8.50) 

The gluon-quark-antiquark antenna function has only one unresolved configuration: if 
quark and antiquark are collinear. In this case, 

£0(1,3,4) ^ —P qq ^ G {z). (8.51) 

«34 

8.4.2 Four-parton antenna functions 

The NNLO real radiation corrections to H — ► gg yield four different four-parton gluon- 
gluon antenna functions: the gluon-gluon-gluon-gluon antenna contains only a leading 
colour (colour-ordered) term F® (l g , 2 g , 3 g , 4 3 ), while the gluon-quark-antiquark-gluon an- 
tennae contain both leading colour and subleading colour contributions G\(l g , 3 q , 4 q , 2 g ) 
and 64(1^,3^,4^,2^). Finally, there is also a quark-antiquark-quark- antiquark antenna 
function (l q , 2 q , 3 q >, 4v), where the quarks are of different flavour. Angular terms are 
indicated where appropriate and will be discussed in detail in Section 8.5 below. 

Like already observed in the quark-gluon case, the colour-ordered gluon-gluon-gluon- 
gluon antenna function F2(l g , 2 g , 3 g , 4 g ) contains several different antenna configurations: 
due to the cyclic nature of the colour indices of the H — > gggg matrix element [47] , each pair 
of two neighbouring partons can represent the two hard partons forming the antenna. In the 
following, the set (ijkl) represents any of the ordered permutations (1234, 2341, 3412, 4123). 

The double unresolved limits are: 
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1. Double soft and soft-collinear limits: 



F4 (1,2, 3,4) 9 Sijki , 

F 4 °(l, 2, 3, 4) ^-^r Sijk Skli , (8.52) 

ff (1,2,3,4) S l , kji {z) - P gg ^ G (z) , 

Sij 

F°(l,2,3,4) ^J^T° S k . Uj (z) — P 99 ^ G (z) • (8.53) 



2. Triple collinear limits: 



F 4 °(l,2,3,4) lsl H k ° P ijk ^ G (w,x,y) . (8.54) 



3. Double collinear limits: 



F°(l,2,3,4) -L- P gg ^ G {z) P gg ^ G (y) . (8.55) 

In the single unresolved limits we have: 

1. Soft limits: 

F 4 °(l, 2, 3, 4) S ljk F°(i, k, I) . (8.56) 

2. Collinear limits: 

F 4 °(l,2,3,4) ^ — P gg ^ G {z) F$((ij),k,l)+wg. . (8.57) 

Sij 

The leading colour gluon-gluon-quark-antiquark antenna function G 4 (l ff , 3 9 , 4^, 2 9 ) is 
colour-ordered: the gluon (l g ) is emitted between gluon (2 g ) and the quark (3 g ), while 
(again due to the cyclic nature of the colour indices) gluon (2 g ) is emitted between gluon 
(l g ) and the antiquark (4^). The second colour-ordering G 4 (l 9 , 4g, 3 q , 2 g ) is obtained by 
exchanging the quark and antiquark momenta. 

The double unresolved limits of this antenna function are: 

1. Double soft and soft-collinear limits: 

G°(l, 3, 4, 2) 3 «^lT° s 12 (3, 4) , (8.58) 



02(1,3,4,2) 3J4 -^ 9 ^° S 1;243 (*) — P qq ^ G (z) , 

G°(l,2,3,4) 3 « I|4 -££T° 5 2;134 W — P q ^ G (z) • (8.59) 

S34 



- 71 - 



2. Triple collinear limits: 

G°(l, 3, 4, 2) 29 W ? P 2 43-gK x, y) . (8.60) 



3. Double collinear limit: 



G°(l, 3, 4, 2) lsl|2 ^ l|4? _J_ P gg ^ G (z) P q ^ G (y) ■ (8.61) 



For the single unresolved limits we have: 
1. Soft limits: 



G°(l,3,4,2) S 213 G°(2,3,4), 

G° 4 (l, 3, 4, 2) ^ 5 124 G°(l, 3, 4) . (8.62) 



2. Collinear limits: 



G°(l,3,4,2) * — P gg ^ G {z) G°((12),3,4) + ang. , 
G°(l,3,4,2) * — P,^ G (z) F 3 °(l,(34),2) + ang. , 
G°(l,3,4,2) * — P^ Q (z) G°(2,(13),4) , 

s 13 

G°(l, 3, 4, 2) 4 *> 9 — P,^q(*) G°(l, 3, (24)) . (8.63) 

In the subleading colour antenna G®(l g , 3 g , 4^, 2 g ) one of the gluons is radiated between 
the quark-antiquark pair(3,j, 4^), while the other is outside the quark-antiquark system. As 
a consequence, no singular structure involves both gluons at once. 

This antenna function contains only one type of double unresolved configuration, the 
triple collinear limit: 

G2(l,3,4,2) U ^ P 3 u-> G (w,x,y) , 

(52(1,3,4,2) 2 ^ 4? P 32 ^ G (w,x,y). (8.64) 
The single unresolved limits always involve one of the gluons. They read: 

1. Soft limits: 

G°(l,3,4,2)^°S 3 i4 G°(2,3,4), 

G2(l, 3, 4, 2) S 324 G° 3 (l, 3, 4) . (8.65) 

2. Collinear limits: 

G° 4 (l,i,j,2)^^P qg ^ Q (z) G§(2,(il),j), 
G°(M,j,2) *M j- P qg ^ Q (z) G§(l,(i2),j) , (8.66) 
where i and j can both play the role of the quark or the antiquark. 
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The antenna function containing two quark-antiquark pairs of non-identical flavours, 
H®(l q ,2q,3 q ',4:qi) contains only one double unresolved limit, the double single collinear 
configuration: 

F 4 °(l, 2, 3, 4) M2 !^' l|4f ' _JL_ P q ^ G {z) P q ^ G {y) ■ (8.67) 

■S12S34 

The two single collinear limits are: 

tf 4 °(l,2,3, 4) * — P M -^ G (z)G§((12),3,4)+ang., 

# 4 °(1,2,3,4) 3q ^' — P q ^ G (z) G°(l,2,(34)) + ang. . (8.68) 
8.5 Angular terms 

Angular terms manifest themselves in collinear limits of antenna functions and matrix 
elements when a final state gluon splits into a quark-antiquark pair or into two gluons. 
Several examples were listed in Section 8. To obtain local subtraction terms (i.e. subtraction 
terms which approach the full multi-parton matrix element in its unresolved limits before 
any integrations are carried out), it is necessary to take proper account of these angular 
terms. In this section, we illustrate an algorithmic procedure to reconstruct the angular 
terms appearing in the simple collinear limits of four-parton antenna functions. The same 
procedure can be generalised in principle to reconstruct angular terms appearing in the 
single and double unresolved limits of multi-parton matrix elements. First steps in this 
direction were performed in [29]. 

In the collinear limits arising from gluon splitting, the four-parton antenna functions 
do not yield the unpolarised splitting functions (8.2) multiplied by a spin- averaged three- 
parton antenna function. Instead, one finds that the four-parton antenna functions factorise 
into the corresponding spin-dependent tensorial splitting functions and tensorial three- 
parton antenna functions [19,40]. 

In constructing the subtraction terms for two colour-connected unresolved partons, 
Section 2.3.2, we used the difference between four-parton antenna function and products of 
three-parton antenna functions (2.17). The former are intended to subtract all singularities 
in the double unresolved region, while the latter ensure that the whole subtraction term is 
free of singularities in all single unresolved regions. If the four-parton antenna appearing 
in (2.17) has simple unresolved limits where angular terms are present, the product of 
two three-parton antenna functions in the same equation will no longer subtract its simple 
collinear behaviour locally. The left-over terms are however vanishing after integration over 
the antenna phase space, thus not affecting the cancellation of infrared poles. 

To access these angular terms correctly, one has to keep track of the transverse mo- 
mentum components of the collinear partons. The collinear limit of partons p\ and P2 is 
defined [40] as the limit k± — ► of 

//,' .; /"' • /.•" "'' 



z 2P-U 1 

pR = (l- z) P» - kft - -^L , (8.69) 
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with 



S12 = 2pi -P2 = ~- 



k\ 



(8.70) 



z(l-z) ■ 

In these equations the vector P^ (P 2 = 0) denotes the collinear direction, while n M is 
an auxiliary light-like vector, which is necessary to specify the transverse component k± 
(k± ■ P = k± ■ n = 0). In the small k± limit (i.e. neglecting terms that are less singular than 
l/k 2 ^), we find the following factorisation formulae of the four-parton antenna functions: 



A°(l,3,4,2 
52(1,3,4,2 
#2(1,3,4,5 
#2(1,3,4,5 
#2(1,3,4,5 
#4°(1,2,3,4 
G°(l,3,4,2 
^2(1,3,4,2 
#2(1,2,3,4 
#2(1,2,3,4 



^ — P^ G (z) (^V(l,(34),2) , 
S34 aa 

|4j 1 



S34 



P£+q(z) (^(1,(34), 2), 



(#sV(1,(34),5) 

S34 aa 

(# 3 V(1,3,(45)) 

S45 aa 
S34 



P^a(z) p 3 V(l,(34),5), 



li 1 



— P^ G (z) (F 3 V((u),fc,0, 

-CoW (GV((12),3,4), 

f W (# 3 V(1,(34),2), 

^-PS^ai') (GV((12),3,4), 

-j-P£^{z) (G°W((34),1,2). 
S34 yy 
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The spin-dependent splitting functions -P' 41 ' appearing in these equations were given in 
[19]. The tensorial three-parton antenna functions (X®)^ u can be derived by analogy with 
the scalar three-parton antenna functions of Sections 5-7 from physical squared matrix 
elements. Their tensorial structure is obtained by leaving the polarisation index of the 
gluon associated with the momentum P^ uncontracted. On contraction of the tensorial 
three-parton antenna functions with a physical gluon polarisation average, we recover their 
scalar counterparts. 

It should be noted that the collinear factorisation (8.71) of four-parton antenna func- 
tions in spin-dependent splitting functions and tensorial three-parton antenna functions 
prevents us a priori from using the scalar (spin-independent) antenna functions derived in 
the previous sections to construct local subtraction terms. It seems that tensorial antenna 
functions are required. However this problem can be circumenvented by explicitly isolating 
the angular terms as follows. 

For each four-parton antenna function X\ yielding angular terms in a given simple 
collinear limit (gluon splitting into two partons i and j in the final state) , one considers an 
angular function denoted by This function must fulfil two properties: (1) it yields 
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the correct local behaviour in this particular limit and (2) it integrates to zero over the 
corresponding unresolved phase space. The second requirement is particularly important, 
since it ensures that the analytic integration will not be modified by the presence of these 
angular functions, and that the integrated scalar antennae are sufficient to fully describe 
the pole structure of the antenna subtraction. In principle, a local counterterm is not even 
required for the numerical implementation; it does however allow a point-by-point check 
of the correct numerical behaviour of the subtraction terms, and guarantees considerable 
improvement of the numerical stability. 

We find that the local counterterms are obtained by the replacement 

X 4 °(l,z,i,2) - x2(l,i,j,2)-@ x o(i,j,z,k ± ) , (8.72) 

where the angular function @ x o(i, j, z, k±) is given by, 



' P£"(z,k ± )(Xi)^(l,(ij),2) -±P t ^ m (z)Xl(l,(ij),2) 



„. . • ■ — - ~-r - - „ 



5.73) 



In this equation, X® is the appropriate three-parton antenna function obtained from X® 
when partons i and j build the parent parton (ij). Pj^u^ stands for the spin-dependent 
splitting function while Pij-ttij) stands for the spin averaged splitting function appropriate 
to the limit under consideration. The collinear momentum fraction z can be expressed in 
terms of invariants formed by the momenta appearing in the four-parton antenna phase 
space. Its precise definition is irrelevant, as long as it yields the correct expression in the 
collinear limit. After the replacement (8.72), the resulting four-parton antenna function is 
locally free from singular terms in the single unresolved regions. 

On the other hand, the term 0^o (i, j, z, k± ) integrates to zero if integrated over the 
unresolved phase space, since 



/ 



d0 k» ± kl f(kl) = -— h\ f(kl) , (8.74) 



with 

pun" _i_ ri^P u 
n-P 

being the gluon polarisation sum in the axial gauge. Applied to the spin-dependent splitting 
function Pj^^y relation (8.74) yields the spin-averaged splitting function 'P^^^y 
Contraction with the tensorial antenna function (X® then reproduces the product of the 
spin-averaged splitting function with the scalar antenna function Pij^(ij) (A3). Integration 
of Q x o over the three-parton antenna phase space made from momenta (i), (j) and either 
(1) or (2) yields zero by construction. 

To illustrate the angular replacement (8.72) on a specific example, we consider the 
four-parton antenna function £?2(1 ? , 3 ? /, 4^/, 2g) in the (3 || 4) limit: 

3, 4, 2) -^2(1,3,4,2) - & A o(i,j, z,k ± ) , (8.75) 
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with 

@ A o(i,j,z,k ± ) = J- P q ^ Q (z,k ± ) (A°V(1,(34),2) - P M -^ (z)A°(l,(34),2) , 

(8.76) 

where z is the momentum fraction carried by one of the collinear partons in the angular 
dependent function Qj^o(i,j,z,k±). 

In the colour-connected double unresolved subtraction term (2.17), B® appears in the 
combination 

B° 4 (l, 3, 4, 2) - E° 3 (l, 3, 4) (34), 2) . 

In the collinear (3 || 4) limit, this expression is not vanishing, but yields some residual 
angular terms. After substituting (8.75), this becomes 



£4(1,3,4,2) 



-j- P%^q{z, k ± ) (A% u (l, (34), 2) - J- P q ^ Q (z)A° 3 (l, (34), 2) 



-£§(1,3,4)A§(1,(34),2), (8.77) 
which is free of singularities in the 3 || 4 limit, since (8.71) holds and, 

£0(1,3,4) 3 -^? -Lp^z). (8.78) 

The term Q A o z, k±) vanishes when integrated over the unresolved phase space, because 
of (8.74). 

9. The 1/N 2 contribution to e + e~ -> 3 jets at NNLO 

To illustrate the application of antenna factorisation on a non-trivial example, in this 
section we derive the l/iV 2 -contribution to the NNLO corrections to e + e~ — > 3 jets. For 
completeness, and also for future reference, we will first discuss all NNLO contributions to 
e + e _ — ► 3 jets in Section 9.1. We construct the 1/N 2 double real radiation subtraction term 
in Section 9.2 and the 1/N 2 virtual single real radiation subtraction term in Section 9.3. 
In Section 9.4, we then show how the integrated subtraction terms cancel the l/iV 2 -poles 
of the two-loop virtual corrections. Some details of the numerical implementation are 
discussed in Section. 9.5. Finally, to illustrate the power of our approach, in Section 10 we 
show that the infrared poles of the two-loop (including one-loop times one-loop) correction 
to 7* — > qqg are cancelled in all colour factors by a combination of integrated three-parton 
and four-parton antenna functions. 

The 0{1/N 2 ) colour or QED-like contribution to three jet production in electron 
positron collisions receives contributions from several different partonic channels. There are 
contributions from 7* — > qqggg and 7* — ► qqqqg at tree-level, 7* — > qqgg and 7* — > qqqq at 
one-loop and 7* — > qqg at two-loops. The four-parton and five-parton final states contain 
infrared singularities, which need to be extracted using the antenna subtraction formalism. 

9.1 The matrix elements 

First we list the tree, one-loop and two-loop amplitudes for 7* — > n partons where n < 5. 
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9.1.1 Tree-level matrix elements for up to five partons 

The tree-level amplitude M®-^ n _ 2 ^ g for a virtual photon to produce a quark-antiquark pair 
and (n — 2)-gluons, 

7* (9) q(pi)q{p2)g{P3) • • • s(?>n) 

can be expressed as sum over the permutations of the colour-ordered amplitude M. A n of 
the possible orderings for the gluon colour indices 

K(n-2) 9 = i<^9) n - 2 E (T a *---T^) lil2 M% n ( Pl ,p 3 ,...,p n , P2 ) . (9.1) 

(i,...,k)eP(3,...,n) 

The squared matrix elements for n = 3, . . . , 5, summed over gluon polarisations, but 
excluding symmetry factors for identical particles, are given by, 



K/ = N 3 A° 3 (l q ,3 g ,2 q ), 



999 
|2 



|2 
99999 1 



iV 4 



E - ^A°(l ? ,3 ff ,4 ff ,2 ? -) 



(ij)€P(3,4) 



(9.2) 
(9.3) 



E ^^A^lqjigjjgjkg^q) Ag ( lg , lg , , kg , 2g ) ^ 



(j,i,fc)GP(3,4,5) 



+ f - /y2 J ^5(!9' 3 9' 4 9' 5 9' 2 9) 



where, 



/°_l 2 

l 9?l 



(9.4) 



(9.5) 



and 



M q / = 4(1 -e)g 2 . 



(9.6) 



The squared colour-ordered matrix elements .A3, and A4 are given in eqs. (5.5), (5.27) 
and (5.28) respectively. For the five parton case, 



A^lg, ig, J 9 , fcg, 2g) |A4^| — 



M° A5 (pi,Pi,Pj,p k ,P2] 



M A5 (pi,Pi,Pj,p k ,P2) + M A5 (p 1 ,p i ,p k ,p j ,p 2 ) + M A5 (pi,p k ,Pi,Pj,P2, 



A5(lq,ig,jg, kg, 2 q ) \M. 



|2 
99 1 



E ■ M ° A5 {pi,Pi,p j , P k,P2) 

(i,j,fc)eP(3,...,5) 



(9.7) 

j 

, (9.8) 
(9.9) 



In the sub- leading colour contribution Ag, gluon is effectively photon-like, while in the 
sub-sub-leading colour contribution (also called Abelian contribution), A°, all three gluons 
are effectively photon-like. Photon-like gluons do not couple to three- and four-gluon 
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vertices, and there are no simple collinear limits as any two photon-like gluons become 
collinear. As a consequence, the only colour connected pair in A® are the quark and 
antiquark. All subtraction terms for this five-parton contribution are therefore based on 
quark-antiquark antennae. 

The tree-level amplitude for 

is given by 

M qqq'q> = * e l5 2 <W<^3<Z4 ^114^3*2 ~ ^ui 2 ^3M^ M% A (pi , p 2 , P3 , Pi) 

+(1^3,2^4), (9.10) 

where S qiq2 S q3q4 indicates the quark flavours. The amplitude M.% 4(pi,P2,P3,P4) thus de- 
notes the contribution from the gi^-pah coupling to the vector boson. The identical quark 
amplitude is obtained 

M° qm = M° qqq , f -(2^A). (9.11) 

The resulting four-quark squared matrix elements, summed over final state quark 
flavours and including symmetry factors are given by 



q,q 
N 4 



N F B° 4 (l q ,3 q Ag,2 q ) - 1 (C 4 °(l 9 ,3 ? ,4 ? -,2 ? -) + C 4 °(2 ? -,4 ? -,3 ? ,l 9 )) 



+N F ^B° 4 (l q ,3 q ,4 q ,2 q ) 
where 

B2(l q ,3 q ,4 q ,2 q )\M° qq \ 2 = \M° B 4 Pl ,P2,p 3 ,P4)\ 2 , 



(9.12) 



C2(l q ,3 q ,4 q ,2 q ) \M° qq \ 2 = -Re(M% A (pi,P2,P3,P4)M ^ 4 (pi,P4,P3,P2^j , (9.13) 
B%(l q ,3 q ,4 q ,2 q ) \M° qq \ 2 = Re[M%4pi,P2,P3,P4)M%(p3,P4,Pi,P2)). (9.14) 



Explicit expressions for B® and C® are given in eqs. (5.37) and (5.42) respectively. The last 
term, B®, is proportional to the charge weighted sum of the quark flavours, Np ;1 , which 
for electromagnetic interactions is given by, 



l^q c q 

It is relevant only for observables where the final state quark charge can be determined. 
There are four colour structures in the tree-level amplitude for 

7*0?) q{pi)q{p2)q' {p?,)q' {P4)g{pb) 
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which reads 

Mqqq'q'g = ^^qiq 2 <*9394 



X 



T iX S isi2 M%° 5 {pi , p 2 , P3 , P4 , P5 ) - T? \ 5i 3 i 4 M°£ 5 (P1,P2,P3,P4,P5) 

1 

AT 

+(1^3,2^4). (9.16) 



^M°B b 5 (pi ,P2,P3,P4,P 5 )~ JjTg u 5 lli2 M^% (pi , p 2 , p 3 , P4 , P 5 ) 



The amplitude A4^ a: 5 (|)i,p2,P3,P4,P5) for x = a, . . . ,d denotes the contribution from the 
(Zi<72~pair coupling to the vector boson. Due to the colour decomposition, the following 
relation holds between the leading and subleading colour amplitudes: 



M°/Jpi,P2,P3,P4,P5) = M°i? 5 (pi,P2,P3,P4, ,P2,P3,P4,P5) 

,P2,P3,P4, ,P2,P3,P4,Pa) • (9.17) 



As before, the identical quark matrix element is obtained by permuting the antiquark 
momenta, 

Km = M° qqqlqlg - (2 <-+ 4). (9.18) 

The squared matrix element, summed over flavours and including symmetry factors is 
given by, 



KJ 2 = EI M ^w/ + Ei M « 

q,q' Q 

= N 5 



N F 

+ 



NN F [B v 5 ' a (l q , 5 g ,4q>;3 q ,,2q) + B^(l q , 4 q/ ;3 q/ ,5 g , 2 q ) 
^ (B° 5 ' c (l q , 5 g ,2 q ;3 ql , 4,-0 + B°/(l q ,2 q ;S q/ , 5 g , 4 q/ ) - 2B° 5 ' e (l q ,2 q ; 3^;5 9 )) 

— C^(l q ,3 q ,4 q ,5g,2 q ) + ^ — ^ ^C^(l q ,3 q ,4 q ,5g,2 q ) + C^(2 q ,4 q ,3 q ,5g, lg^j 

-NN F „ (B° 5 > a (l q , 5 g ,4q>;3 q ,,2 q ) + B°/(l q , 4 q/ ;3 q ,,5 g , 2 q ) - B°' e (l„ 4 q ,- 3 q ,, 2 q , 5,)) 

^> C (1„ 5 g , 2- q - 3 q/ , 4 ql ) + B°/(l q , 2- q - 3 q/ ,5 g , 4- ql ) + B°' c (l„ 2 q ; 3 q ,,4q>\ 5 fl )) 

(9.19) 

where for x = a, . . . , e 

B° 5 ' x (. • •) K/ = \M^ 5 ( Pl , P 2,P3,P4,P 5 )\ 2 , (9.20) 

4°'% • \Kqf = Re (•^Br 5 (Pl^2,P3,P4,P5)M^ 5 ' t (P3,P4,Pl,P2,P5)) , (9-21) 

and 

Cl{l q ,3 q ,4 q ,5 g ,2 q ) |M°/ = -2Re^M ^ 5 ( PuP 2,P3,P4,P5)M ^( Pl ,p 4 ,p 3 ,P2,p 5 ) 
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+ M (P 1 iP2'P3,P4, Pb)M^J (pi,P4,P3,P2,pA 
+ M°z b {pi,P2,PZ,PA,Ps)M°Bj} '(P3,P2,P1,P4,P5) 

+ MB b j5 (pi,P2,P3,P4,P5)M°/j?{p 3 ,P2,Pl,P4,P5)\ 

(9.22) 

C^AqAqAqAg^q) Wqq^ = ~ Re ( M°£ 5 (pi , P2 , P3 , P4 , Ps)M°g£ (pi , P4, P3 , P2 , Pb) ) • 



(9.23) 

9.1.2 One-loop matrix elements for up to four partons 

The renormalised one-loop amplitude M qqg for a virtual photon to produce a quark- 
antiquark pair together with a single gluon, 

7*(g) -> q{pi)q{p2)g{P3) 
contains a single colour structure such that 



Mqqg = ieV2g 



' K.M\ 3 ( Pl ,P3,P2) 



(9.24) 



16vr 2 J ^ li2 " 

Unless stated otherwise, the renormalisation scale is set to [A = q 2 . 

The interference of the one-loop amplitude with the three-parton tree-level amplitude 
(9.1) is given by 



2Re (M°J g M^ g ) = N 3 (J) A^°\l q , 3 g , 2 q ) , 



(9.25) 



where 



4 1X0) (1 9 ,3 9 ,2 5 ) = (n [Al(l q ,3 g ,2g) + Al(s 123 )A° 3 (l q ,3 g ,2g)] 



1 

"N 



Al(l q ,3 g ,2 9 ) + Al(s 123 )A%(l q ,3 g ,2 9 )\ +N F Al(l q ,3 g ,2 q )j , 

(9.26) 

where A\, A\ and A\ are given up to O(e ) in eqs. (5.13), (5.15) and (5.17) respectively. 

The one-loop corrections to 7* — ► 4 partons have been available for some time [8]. The 
one-loop amplitude for 

7* (q) q{pi)q{p2)g{P3)g{p4) 

contains two colour structures, 
Ming = ie2g 2 (yj^) 



E ( Ta ' T<1 % l2 (nm%Api,pupj,P2) - ^mXApi^pj^) 



+\^ a3 5 irt2 M\ d AP^P^P^P2) 



(9.27) 
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where 

M 1 J L d 4 (pi,p 3 ,P4,P2) = M]i d 4 (pi,p4,P3,P2)- (9.28) 

The "squared" matrix element is the interference between the tree-level and one-loop 
amplitudes, 



M°- M l - \ = N 4 ( —) 
qqgg qqgg\ \2-kJ 



(N 2 Al' a (l q ,i g ,j g ,2 q ) -Al'\l q ,i g ,j g ,2 q ) +NN F Al' c (l q ,i g ,j g ,2 q )) 

(ij)eP(3,4) 

_ I (l q ,3g,Ag,2g) — A 4 ' ( 1 g , 3g , 4g , 2g ) ~ — "J A± ( lg , 3g , 4g , 2g ) 



+ -^-^4' 3 g , 4 g , 2 q ) 



N 2 ' 

(9.29) 



where for x = a, . . . , d, 

Al' x (lg,i g ,j g ,2g) \M Q qq \ 2 = Re (m % (pi , Pi , Pj , P 2 )M % (pi ,Pi, Pj ,P2 )) , (9.30) 

Al' x (l q ,3 g ,4 g ,2g) \M°gg\ 2 = Re (^^ 4 (pi,p 3 ,P4,p 2 )>(^ 4 (pi,P3,P4,P2)) , (9.31) 

and 

M\* 4 ( Pl ,p 3 ,P4,p 2 ) = M 1 £ 4 (p 1 ,p 3 ,P4,p 2 ) +M)i* 4 ( Pl ,P4,p 3 , p 2 ). (9.32) 

The renormalised singularity structure of the various contributions can be easily writ- 
ten in terms of the tree-level squared matrix elements multiplied by combinations of the 
infrared singularity operators of Section 4. Explicitly, we find 

Voles(Al' a (l q ,i g ,j g ,2g)) = 2 (l«(e, s u ) + I«(e, Sij ) + lg(e, Sj2 )) A°(l ? , j 9 , 2,), 

(9.33) 

Voles(Al' b (l q ,i g ,j g ,2g)) = 2I q ^(e,s 12 )Al(lg,i g ,j g ,2g), (9.34) 
^e S (4' c (l g , i g ,j g , 2g)) = 2 (l qg ] F (e, su) + I^ F (e, ^) + I« F (e, Sj2 )) i 9 , j„ 2 5 ), 

(9.35) 

Poles(A^ a (l q ,3 g ,4 9 ,2 q )) = 

(21$ (e, S34 ) + Ig(e, ^14) + lg(e, s 23 ) + I$(e, si 3 ) + I$(e, s 24 )) 3 g , 4 g , 2 q ), 

(9.36) 

Vole S (Al' b (l q ,3g,4g,2 q )) = 2I q ^ (e,s 12 )A 4 (l q ,3 g ,4 9 ,2 q ), (9.37) 
Po/e S (I^ c (l„3 3 ,4 3 ,2 g -)) = 

'2I« F (e, S34 ) + 0(6, «i 4 ) + I^ F (e, * 23 ) + I^ F (e, Sl3 ) + lW F (€, S24 )) 

xi°(l„3 ff ,4 ff ,2 ? -), (9.38) 
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roles(Al' d {l q ,3 g A g ,2q)) = 

(21^ (e, s 12 ) + 21$ (e, s 34 ) - 1$ (e,s 14 ) - 1$ (e, s 23 ) - I«(e, Sl3 ) - l£ } (e, S24 )) 

xl°(l„3 ff ,4 ff ,2 9 -). (9.39) 
As at tree-level, the one-loop amplitude for 

7*(g) -> 9(pi)g(p 2 )g'(P3)9 / (P4) 
contains two colour structures, 



( N^X 4 (Pl>P2,P3,P4) - -^A1^ j4 (pi,]72,P3,P4) + N F M £ 4 (pi, p 2 , P3, Pi) 
-T7^li2^3i4 ( ^A1^4(pi,P2,P3,P4) ~ T^^B^CPl: P2, P3, P4) 



+ (1 3, 2 4) , (9.40) 



+iV> X ^4 (Pi > P2 , P3 , PA ) 

where 

Mj? 4 (pi,p 2 ,P3,P4) +M/ A (pi,P2,P3,P4) =M£a(Pi,P2,P3, ,P2,P3,P 4 )- 

(9.41) 

As before, the identical quark matrix element is obtained by permuting the antiquark 
momenta, 

m 9 W = m .W-( 2 ^ 4 )- ( 9 - 42 ) 

Summing over flavours and including symmetry factors, we find that the "squared" 
matrix element, is given by 

2 \M^M\ q \ = 2 \Ki'«' M m'«\ + E 2 K^X 



qqqq 



q,i 



= N4 (— 
V2vr 



NN F Bl' a (l q ,3 q ,4 q ,2 q ) - ^Bl'\l q ,S q ,4 q ,2 q ) + N 2 F Bl' c (l q ,S q ,4 q ,2 q ) 

—C\ ,d (l q , 3 q , 4 q , 2 q ) + C 4 ' 6 (l 9 , 3 q , 4 q , 2 q ) — C\^{\ q , 3 q , 4 q , 2 q ) 
— C4' (2 q , 4q, 3 q , l g ) + ^yf^V ( 2 ?> ^g, lg) Jj-C^ {2 q , 4 q , 3 q , l q ) 

+NN F ^Bl' a (l q , 3 q ,4q, 2 q ) - ^1^/(1,, 3„ 4,-, 2,-) + iV F iV F)7 ^' c (l„ 3„ 4,-, 2,-) 

(9.43) 

where for x = a,b,c 

Bl' x (l q ,3 q ,4 q ,2 q ) \M°qq\ 2 = Re (A^f 4 (pi,P2,P3,P4)A<£J 4 (pi,P2,P3,P4)) , (9.44) 
B 4 1>:e (l„3„4,-,2,-) \M°qq\ 2 = Re (x^ 4 (pi,p 2 ,P3,P4)X^ 4 (P3,P4,Pi,P2)) , (9.45) 
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and for x = d,e, f 

Cl' x (l q ,3 q ,Ag,2g) \M° qq \ 2 = -Re (M^% Pl ,p2,p 3 ,p 4 )M^ A (p 1 ,p 4 ,p 3 ,p 2 )) . (9.46) 

Using the infrared singularity operators of section 4, we can extract the singular con- 
tributions of the renormalised one-loop contribution as, 

Voles(Bl' a (l q ,3 q ,A q ,2g)) = 2 (l$(e, s 14 ) + I$(e, s 23 )) B° 4 (l q , 3„ A q , 2 q ), (9.47) 
Voles(Bl'\l q ,S q ,4 q ,2 q )) = 

2 (2I«(e, Sl4 ) - 21$ (e, s 13 ) + 2I$(e, S23 ) - 21$ (e, s 24 ) + I$(e, a 12 ) + I$(e, a 34 )) 

xB°(l ? ,3 ? ,4 9 -,2 ? -), (9.48) 

PoZes(C]' d (l„3„4,-,2 ff )) = 2 (l$(e, s 13 ) + I$(e, s 24 )) C 4 °(l„ 3„ 4 5 , 2 ? ), (9.49) 

2 («> ai 2 ) + I$(e, su) + iJV (e, s 23 ) + I$(e, s 34 ) - I$(e, s 13 ) - I$(e, s 24 )) 

xC 4 °(l ? ,3 9 ,4 3 -,2 3 -), (9.50) 
Voles(Bl' a (l q ,3 q ,4 q ,2 q )) = 2 (l$(e, s 14 ) + I$(e, s 23 )) 5 4 °(1„ 3„ 4,-, 2,-), (9.51) 

2 (21$ (e, ai4 ) - 21$ (e, s 13 ) + 21$ (e, a 23 ) - 21$ (e, s 24 ) + I$(e, *i 2 ) + I$(e, S34)) 

xB 4 °(l ?1 3 5 ,4 ? -,2,-). (9.52) 

9.1.3 Two-loop matrix elements for three partons 

The renormalised two-loop amplitude M^ g for a virtual photon to produce a quark- 
antiquark pair together with a single gluon, 

j*(q) -> q{pi)q{p2)g{pz) 

contains a single colour structure such that 

/ 2 \ 2 

M^-^^^^J T% 2 M% 3 ( Pl ,p 3 ,p 2 ) . (9.53) 

At NNLO, there are two contributions. One from the interference of the two-loop and 
tree- level amplitudes (9.1), the other from the square of the one- loop amplitudes given in 
(9.24). These terms were computed in [5] and are given by, 

2Re (M^M^) = ^3 (^) 2 Af x0) (l q ,3 g ,2 q ) , (9.54) 

Re (M^XJ = ^3 (^) 2 4 lxl) (l ffJ 3 a ,2 ff ) . (9.55) 

Following [48], we organise the infrared pole structure of the NNLO contributions renor- 
malised in the MS scheme in terms of the tree and renormalised one-loop amplitudes such 
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that, 



Voles (4 2x0) (1, , 3 g , 2g) + 4 lxl) (1„ 3 g , 2 ? )) 

\ ± qqg\ c ) ) £ ± qqg\^> 
,r(l-2e) /A) 



2 



r(l-e) V e 
+2J«( e )4 lx °)(l 9 ,3 g ,2 9 -) 



+ K] jW(2e)+ffg 



(9.56) 



Here, 



J«(e) = iV(4 1 )( e , Sl3 )+4 1 )(6, S23 ) 



■ qg 



qg 



N qq 

+NF(l%, F (e,8 13 )+I% tF (e, a23 j), 



(9.57) 



with the individual defined in (4.1) and 



// 



(2) 

qqg 



0^1 



4eT 1 -e 



4 Cs + -^ 



589 llvr z 



432 72 



, 1, 41 7T 2 \ 

N z + — C3 

* 2 Sd 54 48 / 



+ (- 3C3 "i6 + T M 



19 7T 



18 36 



1 



7T 2 \ iV> 



+ — I NNp + ( - ^ ] + -r-TVf.. 



54 24 J N 27 



(9.58) 



We denote the finite contributions as, 



TiViiW 4( 2x0 Vl q _ Ar2 4 (2x0),/mite .(2x0),/mite 1 ,(2x0),/mite 

j-mtze{j± 3 {L q , Og, z q )) — i\ s± 3jN 2 ^3,1 + 3,i/Af 2 



AT2' 

_l_^Y^^(2xO),/mite Nf_ ^(2x0) , finite 



N 3,Np /N 



+iV 2 4 2x | ) ' / ™ te + iV Fi7 ( ± - AT) Af x J JmUe , (9.59) 



Finite(4 xl \l q ,3 g ,2 q )) = IPA^'"™* + 



2 £ (lxl), finite ^_ A (lx 1), finite _,_ 1 ^(lxl), /mite 



+ 



^2^3,1/^2 



_i_ ?V /\r /t (lxl), /mite ^V"F ,(lxl),/inite 



at2 Alxl), finite 



(9.60) 



Explicit formulae for the finite remainders have been given in [5]. These are expressed 
in terms of one-dimensional and two-dimensional harmonic polylogarithms (HPLs and 
2dHPLs) [50,51], which are generalisations of the well-known Nielsen polylogarithms [52]. 
A numerical implementation, which is required for all practical applications, is available 
for HPLs and 2dHPLs [53]. 
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9.2 Five-parton contribution 

Two different five-parton final states contribute at 1/N 2 to three-jet final states at NNLO: 
7* — > qqggg and 7* — > qqqqg with identical quarks. 

The most colour suppressed contribution to the squared matrix element for 7* — > qqggg 
is proportional to 1/N 2 . It is given by, 

IK-a/ = ^5^g(l„3 fl ,4 s ,5 9 ,2,-), (9.61) 
while the identical quark contribution is 

C§(1„ 2 q , 3„ 4 ? -, 5 fl ) + C 5 °(3 ? , 2,-, 1„ 4 ? -, 5 9 )1 . (9.62) 



|M?J 2 = TV, 



The NNLO radiation term appropriate for the three jet final state is given by 

R 1 j R 

NNLO, A NNLO,C ' 



dcr NNLO - d<T NNLO.A + d(J NNTn C ' ( 9 " 63 ) 



with 



d °NNLO,A = ^d^ 5 (p 1 ,..., P5 -q)^A° 5 (l q ,3 g ,A g ,5 g ,2g)4 5) {p 1 ,...,p 5 ) , (9.64) 
da NNLO,c = fl d< ^' ■ ■ ■ ,P5;q)2Cl(l q ,2 q ,3 q ,4 q ,5 9 )4 5) (pi, • • • ,p 5 ) , (9-65) 

where the symmetry factor in front of A® is due to the inherent indistinguishability of 
gluons. The factor 2 in front of arises from the fact that two different momentum 
arrangements contribute to the squared matrix element (9.62). If the quarks and antiquarks 
are not distinguished by the jet functions, these contribute equally. 

The subtraction terms for these cross sections are straightforwardly obtained using the 
procedure detailed in Section 2. Because there is only a single colour connected pair at 
1/N 2 , there are no singularities associated with almost colour- unconnected or unconnected 
unresolved partons for either of the two contributions. Therefore 

d(T NNLO = da NNLO = °- ( 9 - 66 ) 



The other two types of real radiation subtraction terms are decomposed according to 
(9.63) as 

A ''NNLO — uu NNLO, A ~ r ULI NNLO,C 

The subtraction terms for a single unresolved parton are obtained using (2.16), 



d(J NNLO = da NNT.n A + d(jS Ar r .n n ■ ( 9 - 67 ) 



da NNLo,A = d *s(Pi> ■ ■ ■ >ps; q) I 

i,j,keP c (3,4,5) 
da NNLO,C = ^d*5(pi,...,P5l?)2 

^(l ? ,5 fl ,2 ? -)C°((1J) g ,3 ? ,4 g -,(^)-)4 4) (^ 5 ,p3,P4,^5) 
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+A° 3 (l q ,5 g Aq) C 4 °((15) 9 ,3 ? , (45) 9 -,2 ? -) J^fe^s,^) 
+A°(3„ 5 g , 2 q ) C 4 °(l ? , (35),,% (25)-) J^Cp^Pi^pTs) 
+^(3,, 5 ff ,4 ? -) C 4 °(l ? , (35),, (45) ? ,2 9 -) j| 4) ( P T 5 ,p 1 ,p 2 ,pT 5 ) 
-A° 3 (l q , 5 g , 3 q ) C 4 °((15) g , (35),, 4 ? , 2 ? -) {p&,P2,pt,p£) 

-A° 3 (2 q ,5 g ,4 q ) C2(l„3„ (45),-, (25),-) jffepi,^,^)} • (9.68) 

The sum in the first contribution runs only over the three cyclic permutations of the gluon 
momenta to prevent double counting of identical configurations obtained by interchange of 
j and k. 

The colour-connected double unresolved subtraction terms are given by (see (2.17)), 

da NNLo,A = W5d$ 5 (pi,...,p 5 ;g)^j J2 (A° 4 (i q ,i g ,j g ,2q) 

-A° 3 (l q ,i g ,2 q ) A§((ii),,j ff ,(2i)=) -A° 3 (l q ,j g ,2 q ) AU(M) Q ,i g M)a) 



*Al{J^) q ,k g ,(^)-)jf\p^ j ,p k ^ j ) , (9.69) 
d ^LO,c = ^d<I>5(Pi,...,P5;g)2 

x C2(l„3„4s,2 ff ) A§((134),,5 ff , (234),-) J? \m4,P 5 ,P^) • (9.70) 

In the last equation, no single unresolved terms have to be subtracted from C 4 , which 
contains only a single triple collinear limit, but no single unresolved limits (8.30). 
By construction, the combination 

d&NNLO ~ da NNLO ~ da NNLO (9-71) 

is finite in all unresolved limits. 

The momentum maps used in (9.68)-(9.70) are explicit realizations of the generic 
momentum maps denoted in Section 2. In particular, in the single unresolved subtraction 
term (9.68), we map the three on-shell momenta i, j and k onto the on-shell momenta (ij) 
and (kj) by using the three-to-two-parton momentum map of [25]. With this mapping, 
the single unresolved phase spaces in the collinear limits of the unresolved momentum k 
with the hard momenta i or j are treated symmetrically. This mapping is in contrast to 
the dipole phase space mapping [40], where one collinear limit is mapped out to a larger 
phase space volume than the other. In the double unresolved subtraction terms (9.69) and 
(9.70), we map the four on-shell momenta i, j, k and / onto the on-shell momenta (ijk) 
and (Ijk) by using the ordered four-to-two-parton momentum map of [25]. This mapping 
guarantees that both triple collinear limits of the unresolved momenta j and k with the 
hard momenta i or I are treated symmetrically. The single unresolved limits of j with i 
and of k with I are also mapped in a symmetric manner, ensuring the correctness of the 
double unresolved limit obtained by combining both of them. However, the phase space 
mapping of [25] used here does not yield the correct single unresolved behaviour of the 
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phase space in the simple collinear limits of k with i or of j with I (as explained in detail 
in [25], it is not possible to construct a single four-to-two-parton momentum map which 
yields the correct simple unresolved behaviour in all four limits). We therefore decompose 
A\ into its ordered components (which contain only two of the four simple collinear limits) , 
as explained in the context of (5.28) and use either the ordered mapping onto (ijk) and 
(Ijk) or onto (ikj) and (Ikj). For clarity, these two mappings are not distinguished in 
(9.69). Since C4 contains no simple collinear limit, but only triple collinear limits, it does 
not need to be decomposed any further. 

9.3 Four-parton contribution 

At one-loop, there are two contributions to the colour suppressed contribution proportional 
to 1/-/V 2 , one from the four quark final state and one from the two quark- two gluon final 
state. We therefore decompose the virtual one-loop four-parton contribution as 

^NNLO = d(7 NNLO,A + da NNLO,C ' 

with 

d °NNLO,A = W O d<J, 4(pi,...,P4;g)^if(l g ,3 9 ,4 9 ,2 g -)4 4) ( Pl ,...,p 4 ), (9.73) 

da NNLO,C = j£ (£) d$4(pi,...,P4;g)2C 4 1 ' e (l (? ,3 (? ,4„2 (7 ) Jt\pi,---,P4), (9.74) 

where the origin of the symmetry factors is as in the real radiation five-parton contributions 
of the previous section. 

As detailed in Section 2, there are three subtraction terms; the subtraction of the 
infrared poles in the one-loop contribution, the subtraction of the one-loop unresolved part 
and a compensation term outside the singular region. The first of these three terms is 
constructed in (2.31). Following (2.32), it is minus the integrated form of dcr^ LO given 
in (9.68): 



i VS,l,a N 4 f a s\ Afft i x 1 

d(T NNLO,A = -JV2 [to) d ^(Pl,---»P4;?)2j 

x A%s 12 )A'i(l q ,3 g ,A 9 ,2g)4 4) (pi,P3,P4,P2) , (9.75) 

x {Ai(s 12 )+Ai(su)+A%(s 23 )+Ai(s u )-A%(s 13 )-Ai(s 24 )} 

x C 4 °(l g ,3 (? ,4 9 -,2 g -)4 4) (p 1 ,p3,P4,P2) • (9.76) 

Given that Voles(A 3 (sij)) = — 2/^ (s^), we see that the pole structures of these subtrac- 
tion terms exactly cancel the pole structure of A 4 ' given in (9.37) and of C? e given m 
(9.50). 

Single particle unresolved contributions have to be subtracted only from (9.73), since 
(9.74) vanishes in all its single unresolved limits, implying 

j VS,l,b j VS,l,c n / n 

d(J NNLO,C = d(7 NNLO,C = " ^ 
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Moreover, renormalisation does not affect the colour structure considered here. 
From (2.33), we find: 

, vs,i,b N 4 fa s \ 1 
da NNLo,A = W\2^) d ^(pi,---,P4;g)^ 

•jeP(3,4) V 



+Az(l q ,ig, 2 q ) 



Al((li) q ,j g , (2i) q ) + ^(5l234M§((li),,i fl , (2i) ff ) 

x jf'fe.Pi.pT,) • (9.78) 

The previous term correctly subtracts the singularities in the unresolved limits. How- 
ever away from the limit, where the jet algorithm combines the partons, it introduces 
spurious poles. These terms are compensated by (2.36). In this case there are only two 
hard radiators - the quark and antiquark so that i = n = 1 and j = p = 2: 

VS,l,c N A f®s\ ,^ , N 1 



d(T NNLO,A = Jp {^) d* 4 (pi,...,P4;9)2j 



x AMAl(l q ,i g ,2g)Al((li) q J g ,(2 l )-)jP(n t ,p J ,p^ i ) . (9.79) 

ijeP(3,4) 

S h 

This term corresponds to integrating the last two terms of da^ NLO given in (9.69). 

Taken together, the explicit poles in (9.78) and (9.79) cancel. Explicitly, while dropping 
overall factors of phase space etc., and using the definitions given in Sections 4 and 5, the 
pole part of (9.79) is proportional to 

-21$ (e, s 12 )A° 3 (l q ,i g , 2 q )A°((U) q ,j g , (27)-), 

while the first term in (9.78) is 

(2lg(e, Sl2 ) - 2I$(e,s l2t )) Al(l q J g ,2 q )A° 3 ({U) q ,j g ,(2i) q ). 

The second term in eq. (9.78) is 

24V(€, a i 2i M§(l,,i p ,2,0A§((i0,,j s ,(2i),-)- 

Adding these contributions, we see that the explicit poles cancel. 
By construction, the combination 

, V,l , VS,l,a , VS,l,b , VS,l,c , Q 8n x 

aa NNLO aa NNLO aa NNLO aa NNLO (y.SUJ 

is free of explicit 1/e-poles. Furthermore, it is also finite in its unresolved limits. 
9.4 Three-parton contribution 

The three-parton contribution consists of the two-loop three-parton matrix element to- 
gether with the integrated forms of the five-parton and four-parton subtraction terms, 

, v,2 , s I i VS,1 

aa NNLO aa NNLO + aa NNLO ■ 
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We have shown in the previous subsections that parts of the double real radiation subtrac- 
tion term cancel already with the virtual subtraction term. To make these cancellations 
explicit, we introduce a function daJf NLO , which subtracts these terms from da NNLO and 

V S 1 

adds them to da NNLQ . After these cancellations, we observe that the three-parton contri- 
bution reduces to 



dcr 



V,2 

NNLO 



+ (dc 



NNLO 



dcr 



T 

NNLO, 



+ ( d(J NNLO + & a NNLO 



Upon integration of the unresolved phase space, we find 

1 



d^NNLO ~ & a NNLO 



1 



A%s 12 ) + 2C° 4 (s 12 ) A° 3 (l q ,3 g ,2 q )da 3 , (9.81) 



where we defined the three-parton normalisation factor 

2 



—j d$> 3 (p 1 ,p2,P3;q) J3 (pi,P2,P3) ■ 



(9.82) 



Using the explicit formulae for A® and C4 presented in Section 5, we find that the 
singular contribution is given by 



Voles(da NNLO - da NNL o 



1 

iV2 



2 (l$(e,*i 2 )) - ^(^12) " ^(^12) 
Similarly, integration of (9.78) yields, 

dcr nnlo + dcr NNLO = (Al(si 2 )A%(l q , 3 g , 2 q ) 

+A° 3 (s 12 ) \A\(l q ,3 g ,2 9 ) + Al(s 123 )A° 3 (l q ,3 g ,2 q ) 



A° 3 (l q ,3 g ,2 q )da 3 . (9.83) 



da 3 , 



(9.84) 



resulting in 



Poles(da N S / LO + da NNLO ) = 
4 ( - H^Je, s 12 ) + 21$ (e, s 12 ) [i^l,, 3 9 , 2 ? -) + ^(s 123 )A°(l„ 3 a , 2 5 ) 



da 3 , 
(9.85) 



where it is noteworthy that the finite part of A 3 in (9.84) is cancelled exactly by the terms 
proportional to A\ in A 3 (5.23). 

At 0(1/N 2 ), the singular two-loop virtual contribution is given by expanding (9.56). 
We find, 

Voles(da V / NLO ) = ^2 ( [-2 (/^(e,^)) 2 + 2H^ 1/N2 (e, s 12 )] A° 3 (l q ,3 g ,2 q ) 



+2lg(e, a 12 ) 3 g , 2 q ) + ^( Sl23 )A°(l„ 3 9 , 2 q ) 



da 3 . 



(9.86) 
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Using relation between the hard radiation operators H given in (4.9), it is easy to see 

that 



Voles{da s NNLO + da^NLO + dcr NNLo) = °- ( 9 - 87 ) 

The three-parton contribution is thus free of explicit infrared poles. The remaining finite 
contributions are fully differential in the three-jet variables and can be integrated numeri- 
cally for any given infrared-safe jet definition. 

9.5 Numerical implementation 

Starting from the program EERAD2 [10], which computes four-jet production at NLO, we 
implemented the NNLO antenna subtraction method for the 1/N 2 colour factor contri- 
bution to e + e~ — ► 3j. EERAD2 already contains the five-parton and four-parton matrix 
elements relevant here, as well as the NLO-type subtraction terms da S ^ a NLO and da^^io- 
The implementation contains three channels, classified by their partonic multiplicity: 

• in the five-parton channel, we integrate 

da NNLO ~ dcr NNLO ■ ( 9 - 88 ) 



• in the four-parton channel, we integrate 

da NNLO ~ da NNLO ■ ( 9 - 89 ) 



• in the three-parton channel, we integrate 

da NNLO + d(T NNLO + dcr NNLO ■ ( 9 - 90 ) 



The numerical integration over these channels is carried out by Monte Carlo methods using 
the VEGAS [54] implementation. 

It was already demonstrated above that the integrands in the four-parton and three- 
parton channel are free of explicit infrared poles. In the five-parton and four-parton 
channel, we tested the proper implementation of the subtraction by generating trajec- 
tories of phase space points approaching a given single or double unresolved limit using 
the RAMBO [55] phase space generator. Along these trajectories, we observe that the an- 
tenna subtraction terms converge locally towards the physical matrix elements (it has to 
be pointed out that the 1/N 2 colour structure is free of angular terms, since it does not 
involve any gluon splitting processes), and that the cancellations among individual contri- 
butions to the subtraction terms take place as expected in Section 2. Moreover, we checked 
the correctness of the subtraction by introducing a lower cut (slicing parameter) on the 
phase space variables, and observing that our results are independent of this cut (provided 
it is chosen small enough). This behaviour indicates that the subtraction terms ensure 
that the contribution of potentially singular regions of the final state phase space does not 
contribute to the numerical integrals, but is accounted for analytically. 
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10. Infrared cancellations in e + e — > 3 jets at NNLO 

To finally illustrate the power of the antenna subtraction method at NNLO, we note that the 
infrared poles of the two-loop (including one-loop times one-loop) correction to 7* — > qqg 
are cancelled in all colour factors by a combination of integrated three-parton and four- 
parton antenna functions. 

The integrated five-parton double real radiation subtraction term for three-jet produc- 
tion after cancellation of terms with the virtual single unresolved subtraction term reads: 



with 



& a NNLO ~ & a NNLO ~ X qqg,NNLO ^{l-q, 3 <?> 2 ?) dcJ 3 (l°.l) 



^■qqg^NLO 



N 2 



+N° 



+ N 2 



-\ (AM ~ A° 3 (S12) A° 3 (S 12 )) 

-A 4 (s 12 )-^Al(s 12 )-2C2(s 12 ) 

~AUsu) (V° 3 (s 13 ) + V%s 23 )) + A 3 (s 12 )A%s 12 ) 



^2(*i 2 ) + 2c2(*i2 



+NN F £°(s 13 ) + £%(s 23 ) - X - (Vl(s 13 )£l(s 13 )+V° 3 (s 23 )£l(s 23 )) 



+ 



N F 



Integration of the 
produces 



4 v 

+ \ {V 3 (s 13 )£$(s 23 ) + V$(s 23 )£$(s 13 )) 

- ^(sis) -\S2(S23) - \a° 3 (s 12 ) (£ 3 (s 13 )+£l(s 23 )) 
remainder of the four-parton single real radiation subtraction term 



-Bl{s 12 )- l -£l{ 



.(10.2) 
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with 
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= N 2 
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(10.3) 
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+N F 
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(10.6) 



Adding these integrated subtraction terms to the two-loop three-parton contribution 
to 7* — ► qqg [5], discussed in Section 9.1.3, one observes 



Potes ( dcr^jvio + da^ l LO + dov'tm ) = 



V,2 

T NNLO 



(10.7) 



for all colour structures (9.56). This highly non-trivial cancellation clearly illustrates that 
the antenna functions derived in this paper correctly approximate QCD matrix elements 
in all infrared singular limits at NNLO. The formulae in this section also outline the 
structure of infrared cancellations in e + e~ —> 3j at NNLO, and indicate the structure of 
the subtraction terms in all colour factors. 

The above formulae also illustrate several aspects of the antenna subtraction method 
which were discussed in more generality in the preceeding sections. 

The subtraction terms for two colour-unconnected and almost colour-unconnected un- 

S C S d 

resolved partons da^ NLO (2.24) and da^ NLO (2.26) cancel partially with the virtual one- 
loop single unresolved subtraction term da^^Q (2.36). The pattern of these cancellations 
is different for each colour factor, as can be seen from the different products of two tree-level 
three-parton antenna functions in (10.2). 

In Sections 6 and 8.3.2, we indicated that the tree-level four-parton quark-gluon-gluon- 
gluon antenna function D® contains an unwanted double unresolved limit due to the cyclic 
nature of its colour indices. In the iV 2 -contribution to (10.2), we see how this limit can be 
removed by subtraction of the subleading colour quark-antiquark antenna function A\ and 
its strongly ordered limits. We observe that a similar subtraction has to be carried out on 
the one-loop three-parton antenna function D\ in (10.4). 

11. Conclusions 



In this paper, we presented a new method for the subtraction of infrared singularities in 
the calculation of jet observables at NNLO. We introduced subtraction terms for double 
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real radiation at tree level and single real radiation at one loop based on antenna func- 
tions. These antenna functions describe the colour-ordered radiation of unresolved partons 
between a pair of hard (radiator) partons. 

All antenna functions at NLO and NNLO can be derived from physical matrix elements: 
quark-antiquark antenna functions from the corrections to the process 7* — > qq [44] , quark- 
gluon antenna functions from X ^ 99 [46] and gluon-gluon antenna functions from H — > 
gg [47]. We listed all three-parton tree-level and one- loop antenna functions as well as all 
four-parton tree-level antenna functions in their explicit, unintegrated form. We provided 
all integrals of these functions over the three- and four-parton antenna phase spaces for 
the kinematical situation of a neutral massive particle decaying into massless partons, as 
is relevant to e + e~ — ► m jets. Having the subtraction terms in both unintegrated and 
integrated form, they can now be directly implemented into a parton level NNLO event 
generator for jet production in electron-positron annihilation. 

To demonstrate the application of our new method on a non-trivial example, we im- 
plemented the NNLO corrections to the subleading colour contribution to e + e~ — > 3 jets. 
Starting from an existing programme for e + e~ — ► 4 jets at NLO [10], we constructed the 
newly required subtraction terms for double unresolved real radiation in the five-parton 
channel and for single unresolved real radiation off one-loop matrix elements in the four- 
parton channel. We observed that after implementation of these subtraction terms, five- 
parton and four-parton channels integrate to numerically finite contributions to three-jet 
final states, and that the four-parton channel is free from explicit infrared singularities from 
the loop integrations and the subtraction terms. All explicit infrared singularities obtained 
from integrating the newly introduced subtraction terms cancel in the three-parton channel 
after adding with the the genuine two-loop virtual corrections to this channel, which were 
computed previously [5]. A further illustration of the potential of our method is given in 
Section 10, where we demonstrate that the infrared poles of the two-loop virtual correc- 
tions to e + e~ —>■ 3 jets cancel in all its seven colour factors against linear combinations of 
integrated NNLO antenna functions. 

An immediate application of the method presented here is the calculation of the full 
NNLO corrections to e + e~ — > 3 jets [56]. The antenna subtraction method can be further 
generalised to NNLO corrections to jet production in lepton-hadron or hadron-hadron 
collisions. In these kinematical situations, the subtraction terms are constructed using the 
same antenna functions, but in different phase space configurations: instead of the 1 — ► n 
decay kinematics considered here, 2 — > n scattering kinematics are required, which can 
also contain singular configurations due to single or double initial state radiation. These 
require new sets of integrated antenna functions. 

Acknowledgements 

We would like to thank Gudrun Heinrich and Zoltan Trocsanyi for useful discussions on 
several of the issues presented in this paper and for very constructive criticism expressed 
on our earlier works. 



- 93 - 



We would equally like to thank Alejandro Daleo and Stefan Weinzierl [Al] for pointing 
out a number of typographical errors in the original form of (5.29), (5.30), (5.38), (5.43), 
(6.29), (6.48), (7.13), (7.50), (7.51). These typographical errors are corrected in the present 
version. 

Part of the work presented here was performed during the 2004 "QCD and Collider 
Physics" -programme at the Kavli Institute for Theoretical Physics (KITP), Santa Barbara. 
The authors would like to thank the KITP for its kind hospitality. 

This research was supported in part by the Swiss National Science Foundation (SNF) 
under contracts PMPD2-106101 and 200021-101874, by the UK Particle Physics and As- 
tronomy Research Council and by the National Science Foundation under Grant No. 
PHY99-07949. 

References 

[1] R.K. Ellis, W.J. Stirling and B.R. Webber, QCD and Collider Physics, Cambridge University 
Press (Cambridge, 1996); 

G. Dissertori, I.G. Knowles and M. Schmclling, Quantum Chromodynamics: High Energy 
Experiments and Theory, Oxford University Press (Oxford, 2003). 

[2] E.W.N. Glover, Nucl. Phys. Proc. Suppl. 116 (2003) 3 [hep-ph/0211412]. 

[3] S. Moch, J. A.M. Vcrmaseren and A. Vogt, Nucl. Phys. B 688 (2004) 101 [hep-ph/0403192]; 
A. Vogt, S. Moch and J. A.M. Vermaseren, Nucl. Phys. B 691 (2004) 129 [fiep-ph/0404111]. 

[4] Z. Bern, L.J. Dixon and A. Ghinculov, Phys. Rev. D 63 (2001) 053007 [hep-ph/00 10075]; 
C. Anastasiou, E.W.N. Glover, C. Oleari and M.E. Tejeda-Yeomans, Nucl. Phys. 
B 601 (2001) 318 [hcp-ph/0010212]; 601 (2001) 347 [hcp-ph/0011094]; 605 (2001) 486 
[hep-ph/0101304]; 

E.W.N. Glover, C. Oleari and M.E. Tejeda-Yeomans, Nucl. Phys. 605 (2001) 467 
[hcp-ph/0102201]; 

C. Anastasiou, E.W.N. Glover and M.E. Tejeda-Yeomans, Nucl. Phys. B 629 (2002) 255 
[hcp-ph/0201274]; 

E.W.N. Glover and M.E. Tejeda-Yeomans, JHEP 0306 (2003) 033 [hcp-ph/0304169]; 
E.W.N. Glover, JHEP 0404 (2004) 021 [hep-ph/0401119]; 

Z. Bern, A. De Freitas and L.J. Dixon, JHEP 0109 (2001) 037 [hep-ph/0109078]; JHEP 

0203 (2002) 018 [hep-ph/0201161]; JHEP 0306 (2003) 028 [hep-ph/0304168]; 

A. De Freitas and Z. Bern, JHEP 0409 (2004) 039 [hep-ph/0409007]; 

Z. Bern, A. De Freitas, L.J. Dixon, A. Ghinculov and H.L. Wong, JHEP 0111 (2001) 031 

[hep-ph/0109079]; 

T. Binoth, E.W.N. Glover, P. Marquard and J.J. van dcr Bij, JHEP 0205 (2002) 060 
[hcp-ph/0202266]; 

T. Gchrmann and E. Remiddi, Nucl. Phys. B 640 (2002) 379 [hcp-ph/0207020]; 

S. Moch, P. Uwcr and S. Weinzierl, Phys. Rev. D 66 (2002) 114001 [hep-ph/0207043]. 

[5] L.W. Garland, T. Gchrmann, E.W.N. Glover, A. Koukoutsakis and E. Remiddi, Nucl. Phys. 
B 627 (2002) 107 [hcp-ph/0112081] and 642 (2002) 227 [hcp-ph/0206067]. 

[6] F.V. Tkachov, Phys. Lett. 100B (1981) 65; 

K.G. Chctyrkin and F.V. Tkachov, Nucl. Phys. B 192 (1981) 159; 
V.A. Smirnov, Phys. Lett. B 460 (1999) 397 [hep-ph/9905323]; 



- 94 - 



J.B. Tausk, Phys. Lett. B 469 (1999) 225 [hep-ph/9909506]; 

T. Binoth and G. Heinrich, Nucl. Phys. B 585 (2000) 741 [hcp-ph/0004013]; 

S. Laporta, Int. J. Mod. Phys. A 15 (2000) 5087 [hcp-ph/0102033]; 

T. Gehrmann and E. Remiddi, Nucl. Phys. B 580 (2000) 485 [hcp-ph/9912329]; 

S. Moch, P. Uwcr and S. Wcinzicrl, J. Math. Phys. 43 (2002) 3363 [hcp-ph/0110083]; 

S. Weinzierl, Comput. Phys. Commun. 145 (2002) 357 [math-ph/0201011]; 

C. Anastasiou and K. Melnikov, Nucl. Phys. B 646 (2002) 220 [hep-ph/0207004]; 
V.A. Smirnov, Evaluating Feynman Integrals, Springer Tracts of Modern Physics 
(Heidelberg, 2004). 

[7] Z. Bern, L.J. Dixon and D.A. Kosowcr, Phys. Rev. Lett. 70 (1993) 2677 [hep-ph/9302280]; 
Z. Kunszt, A. Signer and Z. Trocsanyi, Phys. Lett. B 336 (1994) 529 [hep-ph/9405386]; 
Z. Bern, L.J. Dixon and D.A. Kosowcr, Nucl. Phys. B 437 (1995) 259 [hep-ph/9409393]; 

[8] E.W.N. Glover and D.J. Miller, Phys. Lett. B 396 (1997) 257 [hep-ph/9609474]; 
Z. Bern, L.J. Dixon, D.A. Kosowcr and S. Weinzierl, Nucl. Phys. B 489 (1997) 3 
[hep-ph/9610370]; 

J.M. Campbell, E.W.N. Glover and D.J. Miller, Phys. Lett. B 409 (1997) 503 
[hcp-ph/9706297]; 

Z. Bern, L.J. Dixon and D.A. Kosower, Nucl. Phys. B 513 (1998) 3 [hcp-ph/9708239]. 

[9] L.J. Dixon and A. Signer, Phys. Rev. Lett. 78 (1997) 811 [hcp-ph/9609460]; Phys. Rev. D 56 
(1997) 4031 [hcp-ph/9706285]; 

Z. Nagy and Z. Trocsanyi, Phys. Rev. Lett. 79 (1997) 3604 [hep-ph/9707309]; 

S. Weinzierl and D.A. Kosower, Phys. Rev. D 60 (1999) 054028 [hcp-ph/9901277]; 

W.B. Kilgore and W.T. Giele, Phys. Rev. D 55 (1997) 7183 [hep-ph/9610433]; 

Z. Nagy and Z. Trocsanyi, Phys. Rev. Lett. 87 (2001) 082001 [hcp-ph/0104315]; 

Z. Nagy, Phys. Rev. Lett. 88 (2002) 122003 [hcp-ph/01 10315]; Phys. Rev. D 68 (2003) 

094002 [hcp-ph/0307268]; 

J. Campbell and R.K. Ellis, Phys. Rev. D 65 (2002) 113007 [hep-ph/0202176]. 

[10] J. Campbell, M.A. Cullen and E.W.N. Glover, Eur. Phys. J. C 9 (1999) 245 
[hcp-ph/9809429]. 

[11] Z. Bern, L.J. Dixon, D.C. Dunbar and D.A. Kosower, Nucl. Phys. B 425 (1994) 217 
[hcp-ph/9403226]; 

D. A. Kosower, Nucl. Phys. B 552 (1999) 319 [hcp-ph/9901201]; 

D.A. Kosower and P. Uwcr, Nucl. Phys. B 563 (1999) 477 [hcp-ph/9903515]; 

Z. Bern, V. Del Duca and C.R. Schmidt, Phys. Lett. B 445 (1998) 168 [hep-ph/9810409]; 

Z. Bern, V. Del Duca, W.B. Kilgore and C.R. Schmidt, Phys. Rev. D 60 (1999) 116001 

[hcp-ph/9903516]. 

[12] S. Catani and M. Grazzini, Nucl. Phys. B 591 (2000) 435 [hcp-ph/0007142]; 
[13] D.A. Kosower, Phys. Rev. Lett. 91 (2003) 061602 [hcp-ph/0301069]. 
[14] S. Weinzierl, JHEP 0307 (2003) 052 [hep-ph/0306248]. 

[15] C. Anastasiou, Z. Bern, L.J. Dixon and D.A. Kosower, Phys. Rev. Lett. 91, 251602 (2003) 
[hcp-th/0309040]; 

Z. Bern, L.J. Dixon and D.A. Kosower, JHEP 0408 (2004) 012 [hcp-ph/0404293]; 
S.D. Badger and E.W.N. Glover, JHEP 0407 (2004) 040 [hep-ph/0405236]. 



- 95 - 



[16] C. Anastasiou, L.J. Dixon, K. Melnikov and F. Petriello, Phys. Rev. Lett. 91 (2003) 182002 
[hcp-ph/0306192]; Phys. Rev. D 69 (2004) 094008 [hcp-ph/0312266]. 

[17] A. Gchrmann-Dc Ridder and E.W.N. Glover, Nucl. Phys. B 517 (1998) 269 
[hcp-ph/9707224]; 

[18] J. Campbell and E.W.N. Glover, Nucl. Phys. B 527 (1998) 264 [hep-ph/9710255]. 

[19] S. Catani and M. Grazzini, Phys. Lett. B 446 (1999) 143 [hcp-ph/9810389]; Nucl. Phys. B 
570 (2000) 287 [hcp-ph/9908523]. 

[20] F.A. Berends and W.T. Gielc, Nucl. Phys. B 313 (1989) 595; 

V. Del Duca, A. Frizzo and F. Maltoni, Nucl. Phys. B 568 (2000) 211 [hep-ph/9909464]; 
T.G. Birthwright, E.W.N. Glover, V.V. Khoze and P. Marquard, JHEP 0505 (2005) 013 
[hcp-ph/0503063]. 

[21] DA. Kosower and P. Uwcr, Nucl. Phys. B 674 (2003) 365 [hcp-ph/0307031]. 

[22] D. de Florian and M. Grazzini, Nucl. Phys. B 616 (2001) 247 [hep-ph/0108273]; 704 (2005) 
387 [hcp-ph/0407241]. 

[23] A. Daleo, C.A. Garcia Canal and R. Sassot, Nucl. Phys. B 662 (2003) 334 [hep-ph/0303199]; 
A. Daleo, D. de Florian and R. Sassot, Phys. Rev. D 71 (2005) 034013 [hcp-ph/0411212]. 

[24] K.L. Adamson, D. de Florian and A. Signer, Phys. Rev. D 65 (2002) 094041 
[hcp-ph/0202132]; Phys. Rev. D 67 (2003) 034016 [hcp-ph/0211295]. 

[25] DA. Kosower, Phys. Rev. D 67 (2003) 116003 [hcp-ph/0212097]. 

[26] S. Weinzierl, JHEP 0303 (2003) 062 [hep-ph/0302180]. 

[27] W.B. Kilgore, Phys. Rev. D 70 (2004) 031501 [hep-ph/0403128]. 

[28] M. Grazzini and S. Frixionc, JHEP 0506 (2005) 010 [hep-ph/0411399]. 

[29] G. Somogyi, Z. Trocsanyi and V. Del Duca, JHEP 0506 (2005) 024 [hcp-ph/0502226]. 

[30] G. Hcinrich, Nucl. Phys. Proc. Suppl. 116 (2003) 368 [hcp-ph/0211144]; 

C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. D 69 (2004) 076010 
[hcp-ph/0311311]; 

T. Binoth and G. Heinrich, Nucl. Phys. B 693 (2004) 134 [hcp-ph/0402265]; 
G. Heinrich, Nucl. Phys. Proc. Suppl. 135 (2004) 290 [hcp-ph/0406332]. 

[31] A. Gchrmann-Dc Ridder, T. Gchrmann and G. Hcinrich, Nucl. Phys. B 682 (2004) 265 
[hep-ph/0311276]. 

[32] C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. Lett. 93 (2004) 032002 
[hcp-ph/0402280]. 

[33] C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. Lett. 93 (2004) 262002 
[hcp-ph/0409088]; Nucl. Phys. B 724 (2005) 197 [hep-ph/0501130]. 

[34] C. Anastasiou, K. Melnikov and F. Petriello, hcp-ph/0505069. 

[35] R.K. Ellis, D.A. Ross and A.E. Tcrrano, Nucl. Phys. B 178 (1981) 421. 

[36] K. Fabricius, I. Schmitt, G. Kramer and G. Schierholz, Z. Phys. C 11 (1981) 315. 

[37] Z. Kunszt and P. Nason, in Z Physics at LEP 1, CERN Yellow Report 89-08, Vol. 1, p. 373. 



- 96 - 



[38] W.T. Giele and E.W.N. Glover, Phys. Rev. D 46 (1992) 1980; 
Z. Kunszt and D.E. Sopcr, Phys. Rev. D 46 (1992) 192; 

S. Frixionc, Z. Kunszt and A. Signer, Nucl. Phys. B 467, 399 (1996) [hep-ph/9512328]. 

[39] G. Sterman and S. Weinberg, Phys. Rev. Lett. 39 (1977) 1436; 
G. Sterman, Phys. Rev. D 17 (1978) 2773; 17 (1978) 2789. 

[40] S. Catani and M.H. Seymour, Nucl. Phys. B 485 (1997) 291; 510 (1997) 503(E) 
[hcp-ph/9605323]. 

[41] D.A. Kosower, Phys. Rev. D 57 (1998) 5410 [hcp-ph/9710213]; Phys. Rev. D 71 (2005) 
045016 [hep-ph/0311272]. 

[42] F.A. Bcrcnds and W.T. Giele, Nucl. Phys. B 294 (1987) 700; 

M. Mangano, S. Parke and Z. Xu, Nucl. Phys. B 298 (1988) 653; 
M. Mangano, Nucl. Phys. B 309 (1988) 461; 

L.J. Dixon, Proceedings of "Theoretical Advanced Study Institute in Elementary Particle 
Physics (TASI '95): QCD and Beyond", ed. D. Sopcr, World Scientific (Singapore, 1996), 
p.539 [hep-ph/9601359]. 

[43] V. Del Duca, L.J. Dixon and F. Maltoni, Nucl. Phys. B 571 (2000) 51 [hcp-ph/9910563]; 

F. Maltoni, K. Paul, T. Stelzer and S. Willenbrock, Phys. Rev. D 67 (2003) 014026 
[hep-ph/0209271]. 

[44] A. Gchrmann-Dc Riddcr, T. Gchrmann and E.W.N. Glover, Nucl. Phys. B 691 (2004) 195 
[hcp-ph/0403057]. 

[45] A. Gehrmann-De Riddcr, T. Gehrmann and E.W.N. Glover, Nucl. Phys. Proc. Suppl. 135 
(2004) 97 [hep-ph/0407023]. 

[46] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, Phys. Lett. B 612 (2005) 36 
[hep-ph/0501291]. 

[47] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, Phys. Lett. B 612 (2005) 49 
[hcp-ph/0502110]. 

[48] S. Catani, Phys. Lett. B 427 (1998) 161 [hcp-ph/9802439]; 

G. Sterman and M.E. Tejeda-Yeomans, Phys. Lett. B 552 (2003) 48 [hcp-ph/0210130]. 

[49] T. Kinoshita, J. Math. Phys. 3 (1962) 650; 

T.D. Lee and M. Nauenberg, Phys. Rev. 133 (1964) B1549. 

[50] E. Remiddi and J. A.M. Vermaseren, Int. J. Mod. Phys. A 15 (2000) 725 [hep-ph/9905237]. 

[51] T. Gehrmann and E. Remiddi, Nucl. Phys. B 601 (2001) 248 [hcp-ph/0008287]; 601 (2001) 
287 [hcp-ph/0101124]. 

[52] N. Nielsen, Nova Acta Leopoldiana (Halle) 90 (1909) 123; 

K.S. Kolbig, J.A. Mignaco and E. Remiddi, BIT 10 (1970) 38. 

[53] T. Gehrmann and E. Remiddi, Comput. Phys. Commun. 141 (2001) 296 [hcp-ph/0107173]; 
Comput. Phys. Commun. 144 (2002) 200 [hep-ph/01 11255]. 

[54] G.P. Lepage, J. Comput. Phys. 27 (1978) 192. 

[55] R. Kleiss, W.J. Stirling and S.D. Ellis, Comput. Phys. Commun. 40 (1986) 359. 

[56] A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover and G. Heinrich, work in progress. 

[Al] S. Weinzierl, Phys. Rev. D 74 (2006) 014020 [hep-ph/0606008]. 



-97- 



